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ABSTRACT 

Howell et al. have reported the discovery of SN la SN 2003fg (SNLS-03D3bb) 
and conclude that SN 2003fg is very likely a super-Chandrasekhar-mass SN la 
perhaps with a mass ~ 2 Mq. Their work is the first strong evidence that has been 
presented for a super-Chandrasekhar SN la. We have performed an analysis of the 
SN 2003fg data using the Yoon & Langer binding energy formula for a rotating 
super-Chandrasekhar-mass white dwarf (also used by Howell et al.) included 
in a simple model of SNe la (which we call the SSC model for Simple Super- 
Chandrasekhar model for SNe la) which assumes spherically symmetric ejecta 
and relies on the approximations of an exponential density profile for SN la ejecta 
and of a sharp boundary of the SN la iron-peak-element core. Our results support 
the conclusion of Howell et al.: SN 2003fg is very probably super-Chandrasekhar 
and probably has mass ~ 2 Mq. 

Subject headings: supernovae: general — supernovae: SN 2003fg (SNLS-03D3bb) 



1. INTRODUCTION 

Recently, Howell et al. (2006, hereafter H2006), as part of the Supernova Legacy Survey 
(SNLS), have reported the discovery of SN la SN 2003fg (supernova SNLS-03D3bb in the 
SNLS naming scheme). SN 2003fg is a remote supernova a.t z = 0.2440 ± 0.0003 and using 
concordance cosmology {Hq = 70 kms~^ Mpc~^, Qm = 0.3, fiat universe) is found to be 
about 2.2 times brighter in absolute V magnitude than the median brightness of SNe la: it 
is intrinsically the brightest SN la known with any confidence. SN 2003fg's lightcurves in 
shape are not out of the range of normal SN la behavior. Thus, SN 2003fg strongly violates 
the lightcurve width-luminosity relationship (Phillips 1993; Phillips et al. 1999) which is so 
useful in cosmology. 
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The one spectrum published by H2006, which comes from 2 days after B maximum 
(hght) (which H2006 assume to also be bolometric maximum light), looks typical of near-S- 
maximum SN la spectra as H2006's Figure 3 shows by a comparison to the very normal SN la 
SN 1994D (a core normal SN la in the terminology of Branch et al. (2006)). The spectrum 
line features, however, are distinctly narrow. The narrowness suggests somewhat low ejecta 
velocities for the region of line formation in comparison to normal SNe la at the same phase. 
(The width of the lines is determined by the Doppler shift of line opacity in the expanding 
ejecta.) The line velocity of SN 2003fg's Si II A6355 hne (e.g., Wiese et al. 1969, p. 79) (the 
line which is the most characteristic hne of SNe la) in the spectrum is 8000±500 km s~^. (Line 
velocity is the Doppler shift velocity corresponding to the wavelength shift of a P Cygni hne 
absorption trough minimum from the line center wavelength. As supernovae evolve the line 
velocities generally decrease because the ejecta density is falling and the region of sufficient 
opacity to form lines recedes into the ejecta in mass fraction and velocity.) Core normal 
SNe la have photospheric velocities of ~ 11000 km s~^ at 2 days after B maximum (Branch 
et al. 2005, 2006). Since the photospheric velocity is a lower limit on line velocities with 
normal line formation, SN 2003fg's photospheric velocity was unusually low for its phase and 
probably ~ 8000 km s~^ which is, in fact, the photospheric velocity given by H2006 in their 
supplementary information^ 

Note that supernovae after very early times have all mass elements approximately in 
uniform motion, and thus all ejecta structures just scale up linearly with time since explosion 
t. This kind of expansion is called homologous expansion. In homologous expansion, the 
radii r of all mass elements obey 

r ^ vt , (1) 

where v is the mass element velocity and initial radial position in the progenitor is considered 
negligible. Because of equation (1), velocity becomes a good comoving frame coordinate 
for describing supernova ejecta and we conventionally use it as such. Also in homologous 
expansion, all mass element densities scale as t~^. 

The photospheric velocity is the characteristic velocity of the layer of continuum emis- 
sion. In the case of SNe la, the continuum is not a pure blackbody continuum. Because 
of falling density and thus opacity, the photosphere, like the line formation region, recedes 
into the ejecta in velocity coordinate with passing time. As the photosphere recedes, the 
P Cygni lines gradually transform to emission lines. Eventually, the photosphere vanishes 
altogether. The photospheric phase of a supernova is when the photosphere is still important 
to the radiative transfer. After the photospheric phase is the nebular phase. The transition 
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between the two phases is, of course, gradual. 

One point to emphasize is that the SN 2003fg spectrum in the optical was dominated by 
intermediate-mass element (IME) lines and not iron-peak element (IPE) hnes. (We consider 
carbon through calcium as IMEs and scandium through nickel as IPEs.) Optical spectrum 
dominance in the early post-S-maximum phase by IME lines (principally those of Si II, 
S II, and Ca II) is true for all SNe la, not just normal ones. This shows that the spectrum 
formation is not primarily in what we call the IPE core of the ejecta, but at higher velocities. 
Most SNe la show IME-line dominance in the optical spectra in the pre-i?-maximum and 
i?-maximum phases as well. 

The IPE core is the interior region of many SN la explosion models that is almost entirely 
IPEs. At time zero after the explosion the IPE core is in fact mainly radioactive ^^Ni. The 
IPE core is a standard feature of explosion models that result from the standard SN la model. 
In the standard SN la model, the progenitor is a carbon-oxygen white dwarf (a CO WD with 
about equal amounts of C and O) that is very close to the physical Chandrasekhar mass 
(1.38 Mq for CO WDs with equal amounts of C and O: see Appendix A) with typical central 
densities of order 3 x lO^gcm^ and typical central temperatures of order 2.5 x 10^ K (e.g., 
Woosley & Weaver 1994, p. 108). With such densities and temperatures unstable carbon 
burning (i.e., carbon ignition) at or near the WD center will initiate an explosion that totally 
disrupts the WD and creates a SN la. 

The WD is driven to the explosive state by mass accretion from a binary companion. 
The nature of the companion is not certain at present. Two possibilities exist (which are 
not exclusive) for the companion: it could be a post- main- sequence star (or very late phase 
main-sequence star) of some kind or it could be another CO WD that is merging with the 
first WD (e.g., Tornambe & Piersanti 2005, and references therein). The first possibility is 
the single-degenerate (SD) scenario and the second, the double-degenerate (DD) scenario. 
In the SD scenario the accreted matter (which would be mostly hydrogen) burns to CO. 
There are problems with both scenarios. In the SD scenario, the accreted matter may be 
mostly ejected in some way because of unstable helium burning in the accreted matter and 
this prevents the WD central region from reaching explosive conditions (e.g., Piersanti et 
al. 2003, and references therein; Tornambe & Piersanti 2005, and references therein). In the 
DD scenario, it has been found that an off-center carbon ignition may lead to conversion 
of the CO WD to an 0-Nc-Mg WD (e.g., Saio & Nomoto 2004) and perhaps, depending 
on conditions, to collapse to a neutron star rather than a SN la explosion (e.g., Nomoto & 
Kondo 1991). Somehow one of the scenarios or both of the scenarios avoid non-explosive 
fates and allow the CO WD to reach the central conditions that lead to a SN la explosion. 

If a WD gets sufficiently close to the physical Chandrasekhar mass (and does not have 
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high rotation: see below), it should collapse to a neutron star. Electron capture on nuclei 
brings about collapse before the physical Chandrasekhar mass is reached: see Appendix A. 
SNc la happen (at least in some cases) because of carbon ignition in the central regions of 
WDs (as mentioned above) before the point of collapse is reached. 

In a SN la explosion, nuclear burning turns the inner part of the ejecta into nearly 
pure IPEs with the dominant isotope being radioactive ^^Ni: this inner part is the IPE core. 
Outside of the IPE core, IPEs are trace elements. Immediately outside of the IPE core, 
the composition in most explosion models is dominated by explosion-synthesized silicon and 
sulfur. Above that is a layer dominated by explosion-synthesized and pre-existing oxygen 
with probably relatively large abundances of explosion-synthesized silicon and magnesium. 
There may be an outer layer or region of pre-existing CO that survives the burning in some 
SNe la. Carbon lines have been tentatively identified in a few SNe la (e.g., Jeffery et al. 
1992; Fisher et al. 1999; Branch et al. 2003). 

The existence of a somewhat layered composition for SNe la, including in particular 
an IPE core, seems almost certain. Such compositions are the result of partially successful 
near-Chandrasekhar-mass, one-dimensional, hydrodynamic explosion models (e.g., Nomoto 
et al. 1984; Thielemann et al. 1986; Woosley 1991; Woosley & Weaver 1994; Khokhlov et 
al. 1993; Hofiich & Khokhlov 1996; Hoflich et al. 1998) and are indispensable in modeling 
SNe la with synthetic spectra. A reference SN la explosion model of longstanding exhibiting 
a layered composition is the onc-dimcnsional deflagration model W7 (Nomoto et al. 1984; 
Thielemann et al. 1986). In many respects model W7 seems to approximate normal SN la 
behavior to some degree (e.g.. Baron et al. 2006)^. Model W7 does, in fact, have an outer 
layer of unburnt CO (with a trace of Ne) of 0.098 Mq (Thielemann et al. 1986). The unburnt 
mass fraction is 0.071. One-dimensional Delayed-detonation models, which have also been 
partially successful SN la models, typically have much less unburnt mass. For example, 
the near-Chandrasekhar-mass delayed-detonation models of Khokhlov et al. (1993) have an 
unburnt mass fraction of order 0.01 and those of Hoflich et al. (1998) have an unburnt mass 
fraction of order 0.001. Based on this evidence, it is probable that normal SNe la have 
unburnt mass fraction < 0.07. 

A vital part of the standard SN la model is the synthesis of the radioactive ^^Ni in 
the explosion. The ^^Ni is the overwhelmingly dominant source of radiation energy for the 
observable phase of SNe la. The decay chain for ^^Ni is ^^Ni ^^^Co ^^^Fc with half-lives of 
6.077 ± 0.012 days and 77.27 ± 0.03 days for the first and second decays, respectively (e.g.. 



^ Plots of model W7's density profile and composition, and data files that fully specify model W7 are 
available at littp://www.nlin. ou.edu/~jeffery/astro/sne/spectra/model/w7/w7.litml . 
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Firestone & Ekstrom 2004). 

H2006's analysis of the high luminosity of SN 2003fg led to the conclusion that SN 2003fg 
had a very high mass of ^^Ni and a (total) mass of ~ 2M0, and thus that SN 2003fg was 
super-Chandrasekhar (i.e., SNe la with mass greater than the Chandrasekhar mass). H2006 
posit that the progenitor WD had high angular momentum, and thus a high centrifugal force. 
The high centrifugal allowed the progenitor WD to reach a super-Chandrasckhar mass before 
explosive conditions were reached in the central regions. Both SD and DD scenarios may 
allow super-Chandrasekhar SNe la. If the WD is spun up by accretion to very fast differential 
rotation (with mean angular velocity of order a few radians per second on average), then 
the WD may exceed the physical Chandrasekhar mass by up to some tenths of a solar mass 
before reaching explosive conditions in the central region (Yoon & Langer 2005, hereafter 
YL). The simple post-explosion super-Chandrasekhar SN la models discussed in this paper 
(§§ 2 and 4) are independent of whether the pre-explosion evolution followed the SD scenario 
or DD scenario, and so we will not discuss their distinctions further. Before SN 2003fg, there 
was no strong evidence for super-Chandrasekhar SNe la. 

In § 2, we discuss the SN 2003fg ^^Ni mass and (total) mass and H2006's analysis. 
Section 3 discusses the SN la IPE-core boundary and IPE-core velocity which are relevant 
to our analysis of SN 2003fg. Section 4 introduces our SN la model (the SSC model) for the 
analysis of the SN 2003fg ^^Ni mass and mass. In § 5, we compare the predictions of the 
SSC model to observations for a sample of low- 2; SNe la from Benetti et al. (2005) and in § 6, 
to the observations for SN 2003fg. Conclusions and discussion are given in § 7. Appendix A 
discusses the meaning and values of the Chandrasekhar mass. Appendix B presents the 
rotating CO WD binding energy formula of YL which H2006 used and which we use in the 
SSC model and some related results for this binding energy formula. Appendix C shows 
how to solve the SSC model for SN la mass when SN la mass is a dependent parameter. 
Appendix D shows how to solve the SSC model for IPE-core mass when IPE-core mass is a 
dependent parameter. 



2. THE SN 2003fg ^^Ni MASS AND TOTAL MASS 

If we did not know that SN 2003fg was intrinsically exceptionally luminous, it would 
seem like an only somewhat peculiar SNe la and would attract only modest special attention. 
But its luminosity and its violation of the lightcurve width-luminosity relationship makes 
it an exceptionally interesting SN la. In order to produce SN 2003fg's unusual luminosity, 
H2006 have calculated that it had a mass of explosion-synthesized radioactive ^^Ni of Mssni = 
1.29±0.07M«. 
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The calculation for Mssn; made use of the relationship 

M56Ni = -J^ , (2) 

where Lt^oI is the maximum of bolomctric luminosity (which is emitted radiation integrated 
from the ultraviolet to the infrared), £^56Ni(tboi) is the instantaneous rate of release of ra- 
dioactive decay energy per unit mass of ^^Ni by ^^Ni and its daughter isotope ^^Co at the 
time of maximum bolometric luminosity tboi, and a is a correction factor to account for 
time delay in the radioactive decay energy being emitted from the ejecta. H2006's value for 
Lboi/-E'56Ni(tboi) is 1.55 Mq. From analytic solutions for the early SN la lightcurves, Arnett 
(1979, 1982) found that a should be 1 exactly. In realistic SNe la models and in real SNe la, 
one would expect a multiplicity of interacting effects will cause a to be only of order 1 and 
that there will be a range of values with the value for any case depending on that case in per- 
haps a complex way. An early survey of SN la lightcurve calculations suggested a — 1.2 ±0.2 
(Branch 1992), where the uncertainty represents the range of values. An extensive set of 
SN la hghtcurve calculations (Hoflich & Khokhlov 1996) gives a values that extend over the 
range 0.71-1.46 for models that are at least partially successful in fitting lightcurves of one 
or more SNe la: the SNe la included normal SNe la and peculiar SNe la of various sorts. 
The full a range (for models including those that arc not acceptable for known SNe la) is 
0.62-1.46 with mean 1.01 and standard deviation 0.20. There are no unambiguous correla- 
tions of the a values with AfsoNi, rise time to bolomctric maximum light, or model category 
(Hoflich & Khokhlov 1996, Table 2, Fig. 7). However, physically there should be a tendency 
for a to be relatively small if signiflcant ^^Ni matter is near the surface which would happen 
necessarily if the ^^Ni mass was a very large fraction of the ejecta mass. This tendency is 
because some of the radioactive decay gamma-rays can escape the ejecta directly from near 
the surface, and thus not deposit their energy in the ejecta: such escaping gamma-rays are 
not included in the conventional bolometric luminosity. There is a hint of the gamma-ray- 
escape effect in the results of Hoflich & Khokhlov (1996), but it is not clear how to identify it 
among all the other interacting effects. Actually most of the calculated lightcurves of Hoflich 
& Khokhlov (1996) have rise times to bolometric maximum light that are too short for most 
SNe la for which the mean rise time (using B band as a proxy for bolometric luminosity) is 
~ 19 days (Conley et al. 2006). Nevertheless, the calculations of Hoflich & Khokhlov (1996) 
are suggestive of the range of a behavior that may hold for SNe la and we keep that in mind 
for our analysis in §§ 5 and 6. We also will keep in mind though that SN 2003fg is a very 
peculiar and a unique SN la, and so its a value is not well constrained by the a values of 
existing SN la models and other known SNe la. 

Following Branch (1992), H2006 adopted a = 1.2 (which we will consider as the flducial 
value) for all their calculations of Mse^i, but did not consider any variations of a nor include 
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the uncertainty of a in their uncertainty calculations. In fact, their other uncertainties (in 
bolometric luminosity Lboi and rise time to bolometric maximum light ^boi) ^-re negligible in 
comparison to the uncertainty in a. 

Typically, SNe la are thought to have of order 0.6 Mq of ^^Ni. For the sample of well- 
observed 26 \ow-z SNe la studied by Benetti et al. (2005) (hereafter the low- 2; sample), 
H2006 calculate a MseNi range of ~ 0.04-0.86 Mq (using a = 1.2). Thus, H2006's value 
of Mse^i = 1.29 ± 0.07 Mq for SN 2003fg is well outside of the normal Mse^i range and is 
the highest ever calculated with reasonable confidence provided we assume SN 2003fg has 
no extreme variation in a from the fiducial value of 1.2. (We will consider the possible 
variations in a and their effects on (total) SN la mass estimates in §§ 5 and 6.) Note that 
the uncertainties of Lboi/-E'56Ni(tboi) for the low- 2; sample are in some cases greater than for 
SN 2003fg because some of these SNe la are not well out in the Hubble flow where distances 
can be determined with higher accuracy. (Uncertain reddening corrections may also be a 
problem, but reddening corrections were not discussed by H2006.) Thus, the MseNi values 
for the I0W-2; sample are more uncertain than for SN 2003fg. 

The possibihty of SNe la with very high ^^Ni mass has been raised in the past. Assuming 
a long distance for the SN 1991T host galaxy NGC 4527 of 16.4±1.0 Mpc, Fisher et al. (1999) 
suggest a ^''Ni mass of > 1.4 M© for SN 1991T. However, Cepheid distances to NGC 4527 
of 13.0 ± 1.7 Mpc (Gibson & Stetson 2001) and 14.1 ± 0.8 Mpc (Saha et al. 2001) imply 
a less luminous SN 1991T, and thus a lower ^^Ni mass that is in or closer to the normal 
range. For SN 1991T, H2006 find MseNi = 0.86 M© using a = 1.2: it is not clear what 
distance they adopted for NGC 4527. Several SNe la or possible SNe la have very uncertain 
absolute B maxima (Richardson et al. 2002) that are about as bright or even brighter than 
that of SN 2003fg. The most outstanding example is possible SN la SN 19880 for which 
the brightest observed absolute B magnitude is ~ —21.3 (using Hq — 70kms~^ Mpc~^), but 
probably with great uncertainty (Richardson et al. 2002). This B value is ~ 1.2 magnitudes 
brighter than SN 2003fg's absolute B maximum of —20.09 (H2006). The uncertainties are 
such that there is no confidence at present that any of the possibly very bright SNe la are 
as bright or have ^^Ni mass comparable to SN 2003fg. 

The high ^^Ni mass of SN 2003fg has dramatic implications for SN 2003fg's (total) mass. 
As discussed in § 1, in the standard SN la model, a SN la has only about a physical Chan- 
drasekhar mass 1.38 Mq since that is the standard-SN-Ia-model mass of its WD progenitor. 
However, accreting, rotating WDs can exceed the Chandrasekhar mass before reaching the 
density and temperature conditions for explosion. The centrifugal force provides the WD 
with extra support against gravity. YL beheve that WDs of up to ~ 4 Mq are dynamically 
possible, but that possible accretion histories may limit them to not much more than 2 Mq 



-8- 



(Langer et al. 2000). 

Given that the calculated SN 2003fg ^^Ni mass is so close to the Chandrasekhar mass, it 
is difficult on simple grounds to beheve that SN 2003fg's mass is only nearly a Chandrasekhar 
mass. If SN 2003fg were a near-Chandrasekhar-mass SN la and its radiative transfer history 
not exotic, then its photosphere would be well inside its IPE core near B maximum for 
density profiles that work well for normal SNe la. An IPE core of 1.29 Mq of ^^Ni and 
nothing else (which is implausible) in a near-Chandrasekhar-mass SN 2003fg with a model- 
W7 density profile would extend to ~ 15000 km s~^ (Nomoto et al. 1984, Fig. 4). If this 
were the case, the SN 2003fg 2-days-past-i?-maximum spectrum would be dominated by 
IPE lines and, as discussed in § 1, this is not observed. To keep the 2-days-past-i?-maximum 
spectrum looking normal, the actual SN 2003fg IPE-core velocities must probably have been 
kept below ~ 8000 km s"^ (the SN 2003fg photospheric velocity for 2 days past B maximum: 
§ section-intro) by much more than 0.1 Mq of overlying mass. 

A simple scaling argument gives a possible mass for SN 2003fg. Assume that the ratios 
of ^^Ni mass M56Ni and other elements to mass M are roughly constant for SNe la of any 
mass. Then mass will roughly scale with ^^Ni mass. If the average SN la has 0.6 Mq of ^^Ni 
and a near- Chandrasekhar mass, then SN 2003fg would have a mass of 

1 29 



If we now assume that ^J1E|M is approximately a constant as mass is varied, the ejecta 
velocities for the various matter elements of SNe la will be crudely invariant with mass since 
velocities for the various matter elements overall scale approximately as ^J2E/M . Thus, a 
3 Mq SN 2003fg could have the velocities of line formation near B maximum that were not 
wildly different from those of normal SNe la provided the photospheric velocity evolution 
was not wildly different from that of normal SNe la. The assumptions we made to obtain 
the estimate of 3 Mq for SN 2003fg are at most only partially valid, and thus the estimate is 
only suggestive. We verify in Appendix B that ■sj2E/M can be crudely accepted as constant 
with mass M provided super-Chandrasekhar SNe la are approximated by our SSC model 
(see § 4) up to masses of 3Mq. 

Another estimate of the SN 2003fg mass follows from an argument given by H2006 
that we modify a bit here. If an average SN la has a near-Chandrasekhar mass and 
M66Ni = 0.6 Mq, then it has about 0.8 Mq of non-^^Ni matter. If SN 2003fg is to look as 
spectroscopically normal as it does, then perhaps it should have of order the same amount 
of non-^^Ni matter. It follows then that the SN 2003fg mass could be ~ 2.1 Mq. This mass 
estimate is also only suggestive. 

H2006 also estimate the SN 2003fg mass from a simple SN la model that allows for super- 
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Chandrasekhar mass since it assumes a rotating WD progenitor and that has 4 independent 
parameters: (total) mass M, mass of explosion-synthesized IPEs MipE-core; surviving unburnt 
CO mass fraction /co, and progenitor WD central density Pcentrai,WD- H2006 apparently 
in all cases set Pccntrai,WD = 4 x lO^gcm"^ which can be regarded as of order a typical 
progenitor WD central density. We use the symbol MiPE-corc for the explosion-synthesized 
IPE mass since for the SSC model (introduced in § 4), we will make the assumption that 
all the explosion-synthesized IPEs are in the IPE core and no other matter is in the IPE 
core. H2006 do not, in fact, make this assumption although it does not change their model 
much to do so. For the sake of consistent notation, we will use MiPE-core for the description 
of H2006's model. The mass fractions of the ejecta of explosion-synthesized IMEs and 
explosion-synthesized IPEs arc, respectively, fiuE, and /ipe- (Note carbon and oxygen are 
IMEs, but the unburnt amount of CO is not included in /ime-) These are not independent 
parameters since 

/ipe = ^P^core y^^^ _ ^ _ _ j^^^ 

Actually, to obtain model parameter MiPE-core from MseNi, one needs a conversion factor 

^ " M ■ 

^^iPE-core 

H2006 adopt g = 0.7. The rationale for this value is as follows. The deflagration model W7 
gives g = 0.70 (Nomoto et al. 1984; Thielemann et al. 1986) (although to find this number 
accurately one must examine the model W7 data files) and the pure detonation model 
DETl (which is not a reasonable model for any SN la) gives g — 0.66 (Khokhlov et al. 
1993). Models W7 and DETl are rather remote from each other in model space, and so the 
fact that they yield similar g values near 0.7 suggests that g — 0.7 is a good representative 
fiducial value. Obviously g must be in the range to 1. 

Prom equations (2) and (5), the full formula for MipE-core is 

MipE-core = -^^'"[^ , • (6) 

ga±y56^i[tboi) 

As mentioned above, we believe that the uncertainties in the SN 2003fg L^oi and tboi are 
negligible in comparison to the uncertainty of a for SN 2003fg. The uncertainty in g for 
SN 2003fg may be smaller than that of a, but it is still not very well constrained. In our 
analysis in §§ 5 and 6, we will consider how the results change when varying a and g and 
neglect uncertainties in Lboi and tboi- 

The kinetic energy of H2006's SN la model is the energy released by nuclear burning 
minus the binding energy of the WD progenitor and is assumed to be given by 

E — Eip^MlfipE + R{1 — /co — /ipe)] — -E'bind(Pcentral,WD, M) , (7) 
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where burning coefficient E'ipe = LGlfoe/M© (1 foe=10^^ ergs) is tfie energy per unit mass 
of burning CO (witfi equal amounts of carbon and oxygen) to IPEs, M is tlic WD/supernova 
mass in solar mass units, R = 0.768 is the fraction of E'lPE released by burning CO (with 
equal amounts of carbon and oxygen) to IMEs, (1 — fco — /ipe) = /ime from equation (4), 
-E'bind(Pcentrai,WD, M) is the binding energy of the progenitor WD, and Pccntrai.WD is the central 
density of the progenitor WD. The value of ii^ipE is weighted average of the burning coeffi- 
cients for burning CO (with equal amounts of carbon and oxygen) to ^^Ni (l.SGfoe/M©) and 
an average of stable IPEs (1.74 foe/Mo). The weights for the averaging are 0.7 for ^^Ni and 
0.3 for the stable IPEs: these weights agree with the our adoption of fiducial g — 0.7 (see 
above) . The value of R is obtained by dividing an averaged energy per unit mass of burn- 
ing CO (with equal amounts of carbon and oxygen) to IMEs (1.24foe/M0) by l.Glfoc/M©. 
The burning coefficient values 1.56foe/MQ, 1.74foe/MQ, and 1.24foc/Af0 were taken from 
Woosley et al. (2006). (Actually, used -Eipe = l.SSfoe/M© and R = 0.76 which are older 
value given by Branch (1992). The difference between the older and newer values is negligible 
for H2006's analysis.) 

The binding energy formula in equation (7) comes from YL and applies to WDs (in- 
cluding super-Chandrasekhar-mass WDs) partially sustained against collapse by rotation. 
The formula is based on YL's theory of differential rotation of mass-accreting WDs and is 
for zero temperature WDs. The effects of finite temperature are small for the purposes of 
H2006 and the SSC model (which also uses eq. (7): see § 4 and Appendix A). The mean 
angular velocities of super-Chandrasekhar-mass WD models YL consider are typically of 
order a few radians per second (see YL's Tables 1 and 2). The formula is verified for the 
Pcentrai,WD range ~ 10^ — 10^°gcm~^ and the mass range ~ 1.16-2.05 M0. YL's formula and 
rotating super-Chandrasekhar-mass WD models were introduced with the SD scenario as 
the application of primary interest, but they should also apply to the DD scenario if the 
WD in that scenario avoids off-center carbon ignition (Yoon 2006). We present YL's binding 
energy formula and ancillary formulae for reference in Appendix B. 

The Pcentrai,WD range for YL's binding formula includes the range that can lead to WD 
explosion. No explosion probably happens for Pcentrai,WD ^ 2.5 x lO^gcm^^ (e.g., Woosley 
et al. 2006) and above Pcentrai,WD ~ 10^°gcm~^ the WD will probably experience electron- 
capture induced collapse to a neutron star for plausible accretion histories rather than ex- 
plode (Nomoto & Kondo 1991; YL). Woosley et al. (2006) adopt Pcentrai,WD = 2.9 x lO'^ g cm-^ 
as their fiducial explosion density based on their explosion modeling experience: we follow 
them in adopting this density as the fiducial explosion density for most purposes in this 
paper. The i?bind(Pcentrai,WD, -^) IS approximately linear in mass for a fixed Pcentrai,WD (see 
YL's Fig. 9). For each Pccntrai,WD, there are minimum M and -Ebind which are the mass 
and binding energy of a non- rotating WD. The minimum M and £^bind grow slightly with 
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Pcentrai,WD for the range 2.5 x 10^ — 10^°gcm~^: mass from 1.384 M© to 1.414 M©; binding 
energy from 0.501 foe to 0.570 foe (see eq. (B4) and (B2) in Appendix B). Note that the 
non-rotating mass range is above the physical Chandrasekhar mass 1.38 Mq. This is because 
YL use an approximate equation of state for their modehng that gives a Chandrasekhar mass 
of 1.436 (see Appendices A and B). The inconsistency with the physical Chandrasekhar 
mass is small in comparison to the mass variations H2006 and we consider. 

YL's binding energy formula is, in fact, only a close fit to numerical simulations. Those 
simulations rely on many approximations and YL's theory of differential rotation in mass- 
accreting WDs. So YL's binding energy formula is not the last word in binding energy 
formulae. Quantitatively, the binding energy formula is important. If one sets fco = 0.0 
and /iPE = 0.5, which are somewhat realistic for SNe la, then the nuclear burning energy 
part of equation (7) becomes 

Kuc = 1.43 X (M - 1.4) + 2.00 foe , (8) 

where M is solar masses and we have used M — lA Mq as a fiducial value for this argument 
since that is approximately the minimum mass for explosion. The binding energy formula 
can be approximated as a line to high accuracy (Appendix B) with coefficients depending 
on density. For our fiducial Pcentrai.wD = 2.9 x lO^gcm"^, the binding energy formula is 
approximated for the mass range ~ 1.4-3 (which turns out to be the mass range of 
interest in this paper) by the least-squares-fit line 

£;bi„d ~ 1.220 X (M - 1.4) + 0.506 foe , (9) 

where M is solar masses and we have again used M = 1.4 Mq as a fiducial value. Comparing 
equations (8) and (9), it is clear that the binding energy will always be important to the 
behavior of the cjccta for the mass range of interest in this paper (which, again, is ~ 1.4- 
3M0) since the nuclear burning energy is never overwhelmingly dominant for this range. 
Thus, the quantitative results of H2006's SN la model and the SSC model do depend on 
the YL's binding energy formula. But it is likely that improvements to YL's formula will 
not change these results quahtatively, unless rotating WDs with masses significantly greater 
than the Chandrasekhar mass are ruled out by improvements in the input physics of YL's 
theory of differential rotation in mass- accreting WDs. 

Using the kinetic energy from equation (7) with Pcentrai.WD = 4 x lO^gcm"^ in all cases, 
H2006 calculate characteristic kinetic-energy velocities using the formula 
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where the second expression was created using M© = 1.9891 x 10^^ g (e.g., Cox 2000, p. 12). 
H2006's Figure 4b-c shows curves of t'kin versus ^^Ni mass (MseNi) where the curves are 
for WD masses in the range from 1.4 M© (i.e., approximately the Chandrasekhar mass: see 
Appendix A) to 2.1 Mq and a range of fco values from to 0.3. 

For an observation-derived quantity to compare to the V]^in curves, H2006 make use of 
the 26 SNe la of the low- 2; sample. On their Figure 4b-c, H2006 plot Si II A6355 line velocities 
for the low- 2; sample for 40 days past B maximum calculated from a least-squares fit to the 
line velocity evolutions which in some cases extend to 40 days or a bit more past B maximum 
(Benetti et al. 2005). If the composition change between the IPE core and silicon-sulfur layer 
above the IPE core is reasonably sharp (as explosion models suggest: e.g., model W7), then 
the Si II A6355 line velocity of a SNe la should stop decreasing (we can say it plateaus) at 
about the velocity of the IPE-core boundary (which we call the IPE-core velocity for brevity) 
since the line formation for silicon and other IME P Cygni lines cannot recede into the IPE 
core where there is almost no IME matter to cause IME line formation. This plateauing 
of the Si II A6355 line velocity may well have happened in some of the \ow-z sample as 
is suggest by the data from about 40 days past B maximum (Benetti et al. 2005, Fig. 1). 
Unfortunately, line blending with iron lines that develop after B maximum may change the 
Si II A6355 line velocity away from being just the IPE-core velocity. Nevertheless, the long- 
post-5-maximum Si II A6355 line velocities (when the Si II A6355 line can be identified) may 
be a good approximation to the IPE-core velocity. We discuss determination of the IPE-core 
velocities further in § 3. 

H2006 use the slope of the best fit line of that member of the low- 2; sample that is 
closest to SN 2003fg on a plot of near-S-maximum Si II A6355 line velocity versus MseNi 
(H2006's Fig. 4a) to extrapolate from their SN 2003fg 2-day-past-5-maximum Si II A6355 
hue velocity to a 40-days-past-S-maximum Si II A6355 hue velocity for SN 2003fg. They plot 
this on their Figure 4b-c along with the low- 2; sample 40-days-past-S-maximum Si II A6355 
hue velocities. Taking their Figure 4b-c at face value, except for SN 2003fg, all the SNe la 
from the \ow-z sample fall in their figure in the range where a near-Chandrasekhar mass 
is favored or cannot be excluded if /co is high. (The higher /co, the lower /ime low) and 
the less energy from nuclear burning, and so the lower characteristic kinetic-energy velocity 
without having to invoke extra mass beyond the Chandrasekhar mass to slow the ejecta 
down.) If fco is low (or /ime high), some of these SNe la may be super-Chandrasekhar with 
masses up to ~ 2 Mq again taking H2006's Figure 4b-c at face value. 

SN 2003fg is well away from the low-z sample and is clearly in the super-Chandrasekhar 
mass region with mass > 2 Mq for fco ^ 0.2. As discussed in § 1, it is probable that normal 
SNe la have fco ^ 0.07 and, in the context of energy production from the nuclear burning 
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of CO, 0.07 is much less than 0.2. The spectrum of SN 2003fg does show a possible multiplet 
C II A4267 P Cygni line (e.g., Wiese et al. 1966, p. 39) that has hitherto been only very 
tentatively identified in SN la spectra (Branch et al. 2003) and would be expected to be 
weaker in local-thermodynamic equilibrium (LTE) than C II lines that have been tentatively 
identified in SN la spectra. Those other C II lines may be hidden in the noise of the SN 2003fg 
spectra. If the C II line identification is correct, it suggests SN 2003fg may have a larger 
than normal unburnt CO abundance in its outer layers. But given the overall normality of 
the SN 2003fg spectrum, fco ^ 0.2 is probable. Thus, from H2006's Figure 4b-c, one would 
find the SN 2003fg mass to be > 2Mq. Given the conclusion of Langer et al. (2000) that 
WD masses much greater than ~ 2 Mq are not attainable in possible accretion histories, 
H2006's Figure 4b-c suggests a SN 2003fg mass of ~ 2 M©. 

Unfortunately, as H2006 emphasize the comparison of 40-days-past-S-maximum Si II A6355 
line velocities and characteristic kinetic-energy velocity is ambiguous. The characteristic 
kinetic-energy velocity has no definite theoretical relation to the line velocities or the IPE- 
core velocities that the line velocities probably approximate. Thus, the comparison on their 
Figure 4b-c must be regarded as only suggestive. But the fact that SN 2003fg is so far away 
from the other SNe la on the plot and well into the super- Chandrasekhar mass region is 
striking. 

Taken altogether the arguments above give reasonable evidence that SN 2003fg had 
a supcr-Chandrasekhar mass which is perhaps ~ 2Mq. One would like, however, a more 
reliable quantitative value for the SN 2003fg mass or at least more robust lower bounds on 
this mass. In § 4, we introduce a simple model (the SSC model) which is an improvement 
on the model of H2006 and which gives lower bounds on the SN 2003fg mass of varying 
confidence. 



3. THE IPE-CORE BOUNDARY AND THE IPE-CORE VELOCITY 

Before going on to the SSC model, some further discussion of IPE-corc boundary and 
IPE-core velocity (see § 2) needs to be made to support the validity of approximating the 
IPE-core boundary as perfectly sharp in the SSC model (see § 4) and to relate the IPE-core 
velocity to more directly observable quantities. 

The existence of a relatively sharp boundary for the IPE core for normal SNe la at least 
is basically supported by spectrum modeling using explosion models. Typically in useful 
SN la explosion models, the IPEs fall from dominating the composition (nearly 100% of 
the composition by mass) to being a trace (a few percent by mass) over a few thousands of 
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kilometers per second in velocity coordinate. For example, the dominance-to-trace transition 
in the partially successful models W7 (Nomoto et al. 1984; Thielemann et al. 1986) and DD4 
(a well-known delayed-detonation model) (Woosley 1991) happens in the range ~ 8500- 
11000 km s^^. The overall SN la ejecta velocity scale is ~ 25000 km s~^ or more, and so 
assuming a perfectly sharp boundary is reasonable in a simple model. We must add, however, 
that some mixing or explosion-synthesis of IPEs beyond the IPE-cores of explosion models 
is sometimes helpful in SN la spectrum modehng (Baron et al. 2006). We will not consider 
this complication further in this paper it should be kept in mind for future work. 

Since the IPE-corc boundary is not perfectly sharp, one needs a definition for the charac- 
teristic IPE-core velocity. Probably the most sensible definition of the characteristic IPE-core 
velocity for explosion models is the velocity of the spherical shell that would just contain all 
the explosion-synthesized IPEs and nothing else given the model comoving density profile. 
For model W7, this velocity is 9800 km s~^ which just ~ 200kms~^ less than the velocity 
where the mass fractions of ^^Ni and silicon are about equal. For a quick survey, we therefore 
adopt the velocity of equality of ^^Ni and silicon mass fractions as an adequate proxy for the 
IPE-core velocity. For some models that have been found useful in analyzing SNe la, the 
equality velocities fall in the range ~ 7000-13000 km s^^ (e.g., Nomoto et al. 1984; Thiele- 
mann et al. 1986; Woosley 1991; Khokhlov et al. 1993; Hofiich et al. 1998). For the low-;^ 
sample, most of the calculated Si II A6355 line velocities for 40 days past B maximum are in 
the range ~ 8000-1 1000 km s~^ either by measurement or extrapolation: the five faint SNe la 
have lower Si II A6355 line velocities in the range ~ 5000-8000 km s~^ either by measure- 
ment or extrapolation (Benetti et al. 2005). (Faint SN la is a classification of citetbenetti for 
those SNe la that are distinctly fainter than normal SNe la: the prototypes of this classifica- 
tion are SN 1986G and SN 1991bg.) The 40-day-past-B-maximum values for the non-faint 
SNe la arc roughly consistent with equality velocities we cite and support the idea that the 
long-post-B-maximum Si 11 A6355 line velocities at least for non-faint SNe la approximate 
the IPE-core velocities. Recall that some of these 40-day-past-i?-maximum Si 11 A6355 line 
velocities seem to be plateauing (Benetti et al. 2005, Fig. 1). 

Besides the Si 11 A6355 line, other IME P Cygni lines are useful in locating or at least 
constraining the IPE-core velocity. For example, the resonance multiplet Na I A5892 and 
calcium multiplets Ca II A3945 (resonance) and Ca II A8579 (non-resonance, but arising from 
a metastable level) (e.g., Wiese et al. 1969, p. 2, 252, 251) give rise to P Cygni profiles that 
can persist to 100 days or more past B maximum and, in the core normal SN la SN 1994D 
(e.g.. Branch et al. 2005) for example, have line velocities ~ 11000 km s~^, ~ 8000 km s~^, 
and ~ 10000 km s~"^, respectively, at 115 days past B maximum. In the ultraviolet, the 
P Cygni line of resonance multiplet Mg II A2798 (e.g., Wiese et al. 1969, p. 30) also seems 
long persistent (Kirshner et al. 1993; Ruiz-Lapuente et al. 1995) with a line velocity that 



- 15 - 



seems to stay > 14000 km s~^ to 291 days past B maximum although hne-blending with IPE 
hnes makes it hard to be sure that this velocity is representative of the magnesium layer. 
The Mg II A2797.9 region of the ultraviolet is rarely observed for low-^ supernovae (because 
the observation requires space-based telescopes which are not always available), but for high- 
z ones the Mg II A2797.9 region can be redshifted into more readily observed wavelength 
ranges. H2006's spectrum of SN 2003fg almost extends far enough to the blue to see the 
region of the Mg II A2797.9 line. The late, persistent IME P Cygni lines with relatively high 
hne velocities suggest that the line formation layers for these hnes has receded to the vicinity 
of or some point above the IPE core and then plateaued since there is almost no IME matter 
in the IPE core for lower velocity IME line formation. Thus, the sodium, magnesium, and 
calcium line velocities may constrain the IPE-core velocity. On the theoretical side, explosion 
models suggest that the abundances of sodium, magnesium, and calcium all decline going 
inward in the ejecta. For model W7, sodium, magnesium, and calcium fall below their 
solar mass fractions 3.2 x 10~^, 6.0 x 10"'^, and 5.7 x 10^^, respectively, (Asplund et al. 
2005) at velocities ~ 14500 kms^^, ~ 13000 kms^^, and ~ 9000kms~\ respectively. These 
velocities are only representative, of course, but they also suggest that the lowest observed 
hne velocities for Na I A5892, Mg II A2798, Ca II A3945, Ca II A8579 for a SN la could 
constrain the IPE-core velocity. Recall from above, that for model W7 the characteristic 
IPE-core velocity is ~ 9800 km s~^. 

Besides P Cygni hnes, one can also consider for the determination of the IPE-core bound- 
ary the forbidden Fe II and Fe III emission lines that dominate the optical SN la spectrum 
in the nebular phase when the ejecta become optically thin. Modehng these emission hnes 
can be used to constrain the IPE-core velocity. Current results suggest IPE-core velocities 
in the range ~ 7500-10500 km s^^ for normal and somewhat bright SNe la and IPE-core 
velocities of ~ 2000 km s^^ for very faint SNe la like SN 1991bg (Mazzah et al. 1998). Note 
that the calculations for these IPE-corc velocities apparently assumed a uniform density for 
the IPE core (as by Mazzali et al. (1997)) which a significant approximation. 

The above discussion shows both from explosion models and observations that something 
like a sharp IPE-core boundary does exist, a characteristic IPE-core velocity can be defined, 
and that this characteristic IPE-core velocity is in practice determinable to within some 
uncertainty. 



4. THE SSC MODEL FOR SUPER-CHANDRASEKHAR-MASS SNe la 



In this section, we introduce a model for studying super-Chandrasekhar-mass SNe la 
that is an improvement on the model of H2006. We call this model the SSC model for Simple 
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Super-Chandrasekhar model for SNe la. The progenitor for the model is, as in the standard 
SN la model, a mass-accreting CO WD. The SSC model has 3 components. The first 
component is equation (7) (§ 2) (slightly rewritten: see below) to determine the total energy 
of the SN la ejecta. Because we use this equation, we assume we can approximate ejecta 
as consisting of only three types of material for the kinetic energy calculation: explosion- 
synthesized IPEs, explosion-synthesized IMEs, and unburnt CO. By using equation (7), we 
are, of course, also assuming that the progenitor WD is rotating for masses significantly 
greater than the physical Chandrasekhar mass 1.38 Mq (see Appendix A) and we are relying 
on YL's binding energy formula for rotating WDs. As discussed in § 2 our results from the 
SSC model quantitatively depend on YL's binding energy formula. Qualitatively, the results 
from the SSC model are probably robust. 

The second component of the SSC model is the exponential model of the ejecta (Jeffery 
1999) which is explained as follows. Recall from § 1 that after very early times supernovae 
are in homologous expansion where velocity becomes a good comoving frame coordinate 
and all mass element densities scale as t~^. Now as has been known for some time, many 
partially successful SN la models (e.g., model W7) have comoving density profiles that are 
approximately inverse exponentials of velocity, (Jeffery et al. 1992). (Even earlier, it was rec- 
ognized that the comoving density profiles of core-collapse supernova models were sometimes 
approximately piecewise inverse exponentials with velocity (Pizzochero 1990).) Making use 
of homologous expansion and the (inverse) exponential density profile approximation, one 
obtains a useful parameterization of approximate SN la density profiles: 



where to is some fiducial time, pccntrai,o is the central density of the ejecta at the fiducial 
time, and Ve is the e-folding velocity. For examples of SN la model profiles that approximate 
exponentials see, e.g., Dwarkadas & Chevalier (1998, Fig. 1), Baron et al. (2006, Fig. 11), 
and Woosley et al. (2006, Fig. 4). It is straightforward to show that 



where, as in § 2, ii^ is the ejecta kinetic energy and M is the ejecta mass (Jeffery 1999). 
Since SN la kinetic energy should be of order a foe and mass of order Mq, the coefficient 
of the second expression of equation (12) can be take as a good fiducial value for the SN la 
e-folding velocities. To support this we can consider model W7. Using the E and M values 
for model W7 in equation (12) gives Vg — 2670 km s~^ which is not far from 2894.7 km s~^ 
and which, in fact, leads to a good fit to the model W7 density profile: see the URL given in 
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footnote 2 for plots showing this fit. We note that the characteristic kinetic-energy velocity 
defined by equation (10) (§ 2) is related to Vg by 
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It is also straightforward to show that the fraction of ejecta mass, /, interior to reduced 
velocity 
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Note that the exponential density profile extends to infinite velocity, in principal. Ob- 
viously, the approximation of an exponential density profile must fail as velocity becomes 
relativistic. For modehng SNe la, this is not a problem since there is seldom need to invoke 
matter beyond ~ 30000 km s"^ {(5 ^ 0.100) and matter at ~ 40000 km s"^ (/3 ^ 0.133) may 
be the highest ever needed (Jeffery et al. 1992). For such velocities, relativistic effects are 
small since they mostly go as Using our fiducial value = 2894.7km s^^, the reduced 
velocity corresponding to 30000 km s"^ is ~ 10. From equation (15), /(lO) ~ 0.997. Thus, 
in modeling with exponential density profiles, typically much less than 1 % of the ejecta will 
be put at velocities where relativistic effects are important. Therefore, typical relativistic 
effects are going to be negligible when using exponential density profiles. 

As a semi-necessary adjunct to assuming an exponential density profile, we also assume 
that the ejecta can be approximated as spherically symmetric. This is a reasonable sim- 
plifying assumption, but is not necessarily completely valid. Because of their homogeneity, 
normal SNe la are probably quite spherically symmetric: even identical SNe la would look 
significantly different from different viewing directions if they were significantly aspherical. 
There is, however, some evidence from spectropolarimetry for departures from spherical sym- 
metry for both normal SNe la (Wang et al. 2003) and pecuhar SNe la (Howell et al. 2001). 



-18- 



Thus, a priori, one cannot rule out significant asymmetry for any particular SN la especially 
if it is unprecedented like SN 2003fg. Moreover, if the progenitor of SN 2003fg was indeed 
rotating and this accounts for its probable super-Chandrasekhar mass, then its progenitor 
could have been quite aspherical. The ratio of polar to equatorial radii for YL's rotating WD 
models fall as low as ~ 0.3 for masses of ~ 2Mq (see YL's Tables 1 and 2). The explosion 
of such asymmetric WDs may, however, reduce the asymmetry of the matter considerably 
(YL). For the present, assuming spherical symmetry for the ejecta is justified by simplicity 
and the lack of any evidence that SN 2003fg shows any asymmetry. 

The third component of the SSC model is to assume that a well-defined IPE core exists 
with a boundary that we can approximate as absolutely sharp and that we can make the 
approximation that all explosion-synthesized IPE matter is within the IPE core and no other 
elements are. This approximation is reasonable as discussed in § 3. 

Following a similar path to H2006, we take ejecta mass M, IPE-core mass MiPE-core, 
a new parameter /i, and Pcentrai,WD as 4 independent parameters. The new parameter h is 
given by 

, /iME 

^= f ,f ■ (16) 

J CO + JlME 

We prefer h as an independent parameter rather than fco adopted by H2006 or /ime which 
could also be used. One can obtain fco and /ime from h and /ipe (which is obtained from 
equation (4) (§ 2) using M and MipE-core) using the formulae 

/co = (1-/i)(1-/ipe) and /ime = Ml - /ipe) ■ (17) 

Using h as an independent parameter makes simpler the use of IPE-core velocity t'iPE-core as 
an independent parameter in place of M. We want this simplification when we solve for an 
unknown M given MipE-corej 'f^iPE-core) h, and Pcentrai,WD- We find such solutions in § 6 and 
show how the solutions are done and why h preferred to fco and /ime as an independent 
parameter in Appendix C. 

For convenient reference and making use of equations (4) (§ 2) and (17), we rewrite the 
kinetic energy equation (7) (§ 2) in terms of parameter h: 

E = £;ipeM[/ipe + Rh{l - /ipe)] - ^bind(Pcentral,WD, M) . (18) 

Substituting values for M, /ipe, h, and Pcentrai,WD into equation (18), we get the kinetic 
energy E and then substituting values for M and E into equation (12) gives the e-folding 
velocity Ve- 

Now to determine I'iPE-core, we need the reduced IPE-core velocity xiPE-core which we 
obtain by solving for xiPE-core from equation (15) with / set to /ipe- Unfortunately, there is 
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no general analytic solution for x given / for equation (15). In the special cases of / = and 
/ = 1, the solutions are, respectively, a; = and x = oo. Otherwise we can solve equation (15) 
by the Newton-Raphson method (e.g.. Press et al. 1992, p. 355ff) for x given /. The Newton- 
Raphson method is guaranteed to converge in this case since / is monotonically increasing 
with X and only has a stationary point at infinity provided one prevents the Newton-Raphson 
corrections from causing overshooting of the diminishing allowed range for the final solution. 
The derivative of / (which is needed for the Newton-Raphson method solution) is 

df x"^ 

-l = -exp(-x). (19) 

A good initial value for the Newton-Raphson method iteration for / < 1 is the 3rd-order-in- 
small-y solution for x given by 

X3rd = y + ly' + ^y' , (20) 

where y — (6/)^/^. This solution for x is about 10% accurate for / < 0.33: it improves 
in accuracy as y — > 0, of course. Either the 2nd-order-in-small-y or Ist-order-in-small-y 
solutions for x are almost as good initial values for the iteration for any / < 1 as the 3rd- 
order-in-small-y solution for x. As an alternative to the Newton-Raphson method, one can 
solve for x using the iteration equation method using an iteration equation obtained from 
rearranging equation (15). The iteration equation is 



Xout — In 



1 + Xin + (1/2)4 
1-/ 



(21) 



where is an input x value and Xout is an output x value. The iteration is guaranteed 
to converge since the derivative \dxout/dxin\ < 1 with the equality holding only for x^ — 
which is a point that can be avoided since we already know that x = for / = 0. In practice, 
the iteration equation solution to reasonable convergence takes about 50% longer in CPU 
time than the Newton-Raphson method. One can obtain a reasonable approximate analytic 
solution for x, Xapprox, by substituting x^j-d from equation (20) into the iteration equation. 
This approximate analytic solution is given by 



•^approx 111 



l + X3.d + (l/2)xi,d 

1-/ 



(22) 



The approximate analytic solution is always an underestimate with a maximum relative error 
of ~ 19 % for / 0.995: the relative error decreases with / increasing or decreasing from 
/ 0.995 and it goes to zero in the limits / — > and / — > 1. More elaborate analytic 
formulae can give better accuracy for x than equation (22): e.g., substituting Xapprox itself 
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into equation (22) instead of xsrd gives a significantly more accurate approximate analytic 
formula for x. By whatever means, having obtained an accurate xiPE-core; we obtain the 
IPE-core velocity from 

t^IPE-core — ^YPFrcoxe^e ■ 

(23) 

H2006 in their analysis (see § 2) compared 40-days-past-5-maximum-light Si II A6355 
line velocities to Vkin- As we argued in § 2 and § 3, the long-past-S-maximum-light Si II A6355 
line velocity is a reasonable approximation to the IPE-core velocity. As H2006 acknowledged, 
their comparison was ambiguous because there is no necessary theoretical connection between 
the long-past-S-maximum-light Si II A6355 line velocity and t'kin- The situation is further 
clarified by the SSC model, where it is clear that there is no necessity for wipE-core to equal V]^\^ 
even approximately although it could do so. If we set f IPE-core tO fkin (i-*3., ^^jpE-^ore — \/l2), 
then /ipE would be fixed at 0.6724619703381 ~ 2/3 as we see from equation (15). However, 
in the SSC model we are free to vary /ipe from to 1 provided that E >Q which is necessary 
for the WD to explode. 

We illustrate the behavior of the SSC model in Figures 1 and 2 by plotting "(^iPE-corc 
as a function of MipE-corc foi' representative values of the 3 other independent parameters 
(i.e., mass M, /i, and Pcentrai.wo) which are held constant for each curve. (The data points 
on the figures are discussed in §§ 5 and 6.) Recall from § 2 that the Pcentrai,WD range 
~ 2.5 X 10^ — lO^^gcm"^ is believed to be the allowed range for WD explosions and that 
our fiducial Pcentrai,WD = 2.9 X lO^gcm"^. To explore the possible Pcentrai,WD range for SN la 
progenitors, we used the fiducial Pcentrai,WD = 2.9 x lO^gcm"^ (which is not far from the 
low end of the range) for the Figure 1 calculations and Pcentrai,WD = 10^°gcm^^ for the 
Figure 2 calculations. Most explosion modelers find Pccntrai,WD much closer to the low end 
of the allowed range for plausible explosions (e.g., Nomoto et al. 1984; Thielemann ct al. 
1986; Khokhlov et al. 1993; Woosley & Weaver 1994; Hofiich & Khokhlov 1996; Hofiich et 
al. 1998): thus Pcentrai,WD = lO^'^gcm"^ is actually an extreme upper limit and a WD that 
evades an explosion and gets to this limit is probably on the verge of collapse to a neutron 
star (Nomoto & Kondo 1991; YL). The binding energy increases with central density (see 
YL's Fig. 9), and so all other things being equal, increasing density leads to lower fipE-core- 
This is why the curves in Figure 2 are generally lower than their counterparts in Figure 1. 

The solid curves (which form cornucopias) are ior h — 1 which means /ime = 1 — /ipe, 
/co = 0, and there are only IMEs above the IPE core. The mass in solar mass units for each 
solid curve is given at the upper right end of the curve: the masses run from 1.388 Mq for 
Figure 1 and 1.414 for Figure 2 to 3Mq in increments of 0.1 M© for M > 1.5 M©. The 
masses 1.388 Mq and 1.414 Mq are the masses for non-rotating WDs from YL's equation (22) 
(see also eq. (B4) in Appendix B) for the input Pcentrai.WD values and we just label their curves 
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by Mnr on the figures. The non-rotating WD masses are the lower bounds on the masses 
for the input Pccntrai,WD values and are near the physical Chandrasekhar mass 1.38 (see 
Appendix A). As mentioned in § 2, the non-rotating WD masses can slightly exceed the 
physical Chandrasekhar mass 1.38 Mq. This is because YL use an approximate equation of 
state for their modeling that gives a Chandrasekhar mass of 1.436 M0 (see Appendices A 
and B). This inconsistency from the physical Chandrasekhar mass is small in comparison 
to the mass variations we are considering. Actually all the masses we discuss probably have 
uncertainties of a few percent for the given input parameters because they are based on 
YL's choice of equation of state: they are probably mostly shghtly larger than the input 
parameters would imply with a more accurate equation of state (see Appendix A). For 
clarity, we call the 1.388 and 1.414 mass curves the Mnr mass curves. 

The general behavior of the solid curves can be explained as follows. In regard to varying 
mass M for a given IPE-core mass, the smaller M is, the higher t'lPE-core tends to be since the 
WD binding energy decreases with decreasing M, there is less mass to accelerate, viPE-core 
moves to a higher mass fraction in the ejecta, and the nuclear burning energy decreases 
rather weakly only through decreasing amount of burning to IMEs. This explains why as 
curve mass decreases, the curves move upward. In overall curve appearance, the curves 
move toward the upper left as mass decreases. In regard to varying MipE-core ioi fixed M, 
as MiPE-core iucreases, viPE-core increases monotonically. There are two causes for this. First, 
as one moves outward in mass fraction with everything fixed, ejecta velocity will increase 
since outer matter must move faster than inner matter. Second if as one moves outward 
in mass fraction and the IPE-core boundary moves with you, there is more kinetic energy 
than otherwise to drive the outer layers. Both causes help viPE-core to increase monotonically 
with MiPE-corc- Note that the solid curves all begin at AfiPE-corc = since burning all the 
non-IPE-core matter to IMEs always ensures adequate energy for a WD explosion. The solid 
curves also begin at vipE-core = since that is the limiting IPE-core velocity when there is 
no IPE core. 

The dotted curves are for /i = which means /ime = 0, fco = 1 — /ipe, and there is only 
unburnt CO above the IPE core. As mentioned in § 1, the spectra of all SNe la, including 
SN 2003fg, are dominated by IME lines at least in the early-post-5-maximum phase, and 
so we know that there must be abundant IMEs in the outer layers of all SN la ejecta. 
Partially successful explosion models also yield abundant IMEs: for example. Model W7 
gives /ime = 0.281 and h — 0.798 (Thielemann et al. 1986); example delayed-detonation 
models give h > 0.95 (Khokhlov et al. 1993; Hofiich et al. 1998). Thus, must be regarded 
an extreme lower limit on h. However, as discussed in § 2, SN 2003fg may have had more 
unburnt CO than normal SNe la, and so h may be significantly less than 1 for SN 2003fg. 
Thus, the h — curves are relevant to this study. 
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The masses for the dotted curves are the same as for the sohd curves and dotted curves 
are ordered in the same way as the sohd curves by mass: the dotted curves are labeled by 
mass at their low IPE-core mass end. On Figure 2 the dotted curves for masses greater 
than 2.8 Mq arc off the figure. Since there is no burning of IMEs for the dotted curves, they 
represent models with less energy than the counterpart solid curves. Thus, fiPE-corc values 
for the dotted curves at a given MiPE-core are always less than that for the counterpart solid 
curves at the same MipE-core- The dotted curves all begin at zero IPE-core velocity and 
at non-zero IPE-core mass. In the case of the dotted curves, the IPE core must reach a 
sufficient size before there is enough nuclear burning energy to unbind the WD: recall there 
is no burning of matter outside of the IPE core in the case of these curves. The behaviors 
of the dotted curves as M is varied for fixed MiPE-core and as MiPE-core is varied for fixed M 
have essentially the same explanation as for the solid curves. 

As MiPE-core grows for a constant M, the relative amounts of IMEs and CO become less 
important for the energetics, and thus the dotted curves converge toward their counterpart 
solid curves as MiPE-core grows large. For the counterpart dotted and solid curves for larger 
M, close convergence occurs off the top or right edges of the figures. Final convergence, of 
course, happens in the limit of MipE-core M which causes fipE-corc to go to the unphysical 
limit of infinity. Non-relativistic physics fails long before this velocity limit is reached. Since 
the figures only show up to 13000 km s~^ (which \s (3 ^ 0.0434), the region where non- 
relativistic physics fails occurs well off the figures for all the curves. 

Wc need to emphasize that YL have only verified their binding energy formula for the 
mass range ~ 1.4-2.05 M©, and thus we have made a large extrapolation of the usage of 
binding energy formula in plotting curves for masses up to 3M0. 

5. COMPARING THE SSC MODEL AND OBSERVATIONS OF THE 

LOW-z SAMPLE OF SNe la 

The circle points on Figures 1 and 2 are for the low- 2; sample of SNe la from Benetti et 
al. (2005) which was discussed in § 2. The + points are discussed below, and the letters N, 
O, X, and Y are for SN 2003fg and are discussed in § 6. The IPE-core masses for the circle 
points were determined using equation (6) with the L-i,o\/ Ei&m{th6\) values of H2006, a — 1.2 
(as adopted by H2006 and which we take to be fiducial for this paper), and g — 0.7 (which, 
following H2006, we argued was a good representative fiducial value in § 2). The IPE-core 
masses are in the range ~ 0.04-0.86 Mq. The velocities are the Si 11 A6355 line velocities from 
35 days past B maximum derived using the least-squares line fits of Benetti et al. (2005). 
The velocities are in the range ~ 5900-11200 km s~^. As mentioned in § 2, the uncertainties 
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of Lboi/-E'56Ni(iboi) for the low- 2; sample are in general greater than for SN 2003fg because 
some of these SNe la are not well out in the Hubble flow. We will neglect this complication 
in our analysis for simplicity and since we are primarily just interested in overall behavior 
of the low- 2; sample when analyzed using the curves on Figures 1 and 2. 

We prefer 35 days past B maximum to H2006's 40 days past B maximum in calculating 
the Si II A6355 line velocities to be used in comparison to the SSC model curves. The 
reasoning is that after about 35 days past B maximum the derived values in most cases 
would be extrapolations of uncertain quality of the Si 11 A6355 line velocity curves (Bcnctti 
et al. 2005, Fig. 1). Also some of the Si II A6355 line velocity curves seem to be platcauing 
at about 35 days past B maximum (Benetti et al. 2005, Fig. 1) which suggests that the 
Si II A6355 line formation layer has stopped receding into the ejecta because it has receded 
to the boundary of IPE core and there is no significant sihcon below that layer. Thus, 
the 35-days-past-S-maximum values may in many cases best approximate the IPE-core 
velocities and this is what we hope the long-past-S-maximum Si II A6355 line velocities will 
approximate. 

The Si II A6355 line velocity of SN 1991bg may still be decreasing at 35 days past 
B maximum which is about as far as the data of Benetti et al. (2005) go, and so SN 1991bg's 
35-day-past-i?-maximum Si II A6355 line velocity is only an upper limit on the IPE-core 
velocity. From modeling, Mazzali et al. (1998) found that the IPE-core velocity of SN 1991bg 
was ~ 2000 km s^-*^. SN 1991bg is one of the 5 SNe la designated as faint by Benetti et al. 
(2005) in the low- 2; sample. These faint SNe la have the lowest IPE-corc masses, and so give 
the leftmost circle data points on Figures 1 and 2. In order of increasing IPE-core mass on 
Figures 1 and 2, the faint SNe la are: SN 1999by, SN 1991bg, SN 1997cn, SN 1986g, and 
SN 1993H. The situation for the Si II A6355 line velocities of the other faint SNe la are not 
as clear as for SN 1991bg, but their day-35-past-S-maximum Si II A6355 line velocities are 
more likely to be only upper bounds on the IPE-core velocities than is the case for the other 
SNe la in the low- 2; sample. 

The circle points are rather dispersed and do not he along the Mnr mass curves. Some 
are to the right of the Mnr mass curves which suggests the possibihty that they are for 
slightly super-Chandrasekhar SNe la. But this is very uncertain. The circle points to the left 
of the Mnr mass curves are not to be taken for (significantly) sub-Chandrasekhar SNe la. A 
mechanism for sub-Chandrasekhar SNe la has been proposed (Livne & Glasner 1990, 1991), 
but this mechanism has been ruled out. Explosion models produced using the mechanism 
led to synthetic lightcurves (Hoflich & Khokhlov 1996) and spectra (Nugent et al. 1997) that 
failed to match observations: among other things the synthetic lightcurves and spectra were 
too blue. Thus, the circle points to the left of the Mnr mass curves are all likely to be 
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from near-Chandrasekhar (or less likely somewhat super-Chandrasekhar) SNe la. The main 
explanation for the deviations of the circle points from the Mnr mass curves is likely that a 
is not 1.2 in general. The faint group of SNe la may also deviate significantly from the Mnr 
mass curves because their day-35-past--B-maximum Si II A6355 line velocities arc not close 
upper bounds on their IPE-core velocities. Of course, the SSC model may also be somewhat 
inadequate for real SNe la and/or the g factor may vary significantly from the fiducial 0.7 
value. 

The reason for believing it is varying a that leads to the dispersion of the circle points 
from the A^nr mass curves is, as discussed in § 2, that there is reason to believe that a can 
have a large range: recall that the modeling results of Hofiich & Khokhlov (1996) suggest a 
range of 0.62-1.46 for a. We can derive a semi-empirical range for a, by varying a for the 
low-z sample such that all the points lie on the Mnr mass curves. The range is semi-empirical 
because we are relying on observations and on the SSC model in the non-rotating mass limit. 
The range is based on the assumption that most SNe la have nearly the Chandrasckhar mass 
(which is approximated by Mnr in the SSC model). Wc will only do the variation for the 
Mnr mass curves for /i = 1 since this is likely to be closer to the truth for normal SNe la 
than curves for h = (see § 4). The + points on Figures 1 and 2 are the result of varying a. 
Note that a < 1.2 moves points to larger IPE-core mass and a > 1.2 moves points to smaller 
IPE-core mass. The range for a is 0.13-1.72 for Figure 1 and 0.13-1.67 for Figure 2. If we 
exclude the five faint SNe la for which the day-35-past-S-maximum Si II A6355 hne velocities 
are more hkely to be only upper bounds on the IPE-core velocities than for the other SNe la, 
we find the a ranges to be 0.68-1.72 for Figure 1 and 0.65-1.67 for Figure 2. Given the 
assumptions that the SSC model is valid, the non-faint SNe la in the low-^ sample have 
day-35-past-B-maximum Si II A6355 line velocities are nearly the IPE-core velocities, the g 
factor is 0.7, /i ~ 1, and that none of the low- 2; sample is significantly super-Chandrasekhar 
(nor sub-Chandrasckhar), then the semi-empirical a range can be taken as ~ 0.6-1.7. This 
a range is plausible given the results of Hofiich & Khokhlov (1996), but it is not definitive, 
of course. If we used the h — curves instead of the h — 1 curves, the semi-empirical range 
for a would shift to lower values. 



6. COMPARING THE SSC MODEL TO THE SN 2003fg OBSERVATIONS 

We now turn our attention to comparing the SSC model to the SN 2003fg observations. 
In the case of SN 2003fg, we cannot adequately determine the day-35-past-i?-maximum 
Si II A6355 line velocity given only the near-S-maximum Si II A6355 line velocity of 8000 ± 
500kms~^. Therefore, we will only adopt 8000 km s~^ as an upper bound on the IPE-core 
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velocity. The points N, O, X (partially overlapping with O), and Y on Figures 1 and 2 are all 
SN 2003fg IPE-core velocity upper bound points of 8000 km s~^. As Figures 1 and 2 show, 
since we can only set an upper bound on the SN 2003fg IPE-core velocity, we can only set 
a lower bound on the SN 2003fg mass. 

Unfortunately in trying to set such a lower bound on SN 2003fg mass, there is no 
secure way to estimate the IPE-core mass MipE-core (which is a SSC model parameter) 
primarily because of the dispersion in a and secondarily because we can only estimate the 
g factor. Recall Mipg-core is calculated from equation (6) in § 2 using g and a and H2006's 
Lboi/-£'56Ni(tboi) value 1.55 M0 which we assume to have negligible uncertainty The values 
of the SSC model parameters Pccntrai.wD and h are also not fully constrained. Therefore, we 
calculate lower bound masses varying the four parameters Pcentrai,WD, 5', and h from what 
we believe to good representative fiducial values to what we believe to be extreme values 
which tend to minimize the lower bound masses. (Recall, they will be lower bound masses 
because 8000 km s~^ which we use as the input value for vipE-core is only an upper bound 
on f IPE-core-) The fiducial values are Pcentrai,WD = 2.9 x lO^gcm"^ (which we adopted as 
fiducial in § 2 following Woosley et al. (2006)), (7 = 0.7 (which is certainly good for some 
SN la models: see § 2), a = 1.2 (which is good in the sense that it is not extreme: see 
§§ 2 and 5) and h = \ (which is good for normal SNe la given that their spectra up until 
after B maximum are dominated by IME fines and given that model W7 has h — 0.798 
(Thielemann et al. 1986) and some example delayed-detonation models have h > 0.95 (e.g., 
Khokhlov et al. 1993; Hoflich et al. 1998)). The extreme parameter values we consider are 
Pcentrai,WD = 10^° g cm~^ (which is extreme based on the fact the WD is on the verge of 
collapse rather than explosion for this value: Nomoto & Kondo (1991); YL), g = 1 (which 
is an absolute upper bound: sec § 2), a = 1.7 (which is the upper end of the semi-empirical 
range for a that we derived § 5 and is larger than the largest a value of Hoflich & Khokhlov 
(1996) (see § 2)), and h = (which is extreme for normal SNe la given that spectra until after 
5-maximum are dominated by IME lines and for SN 2003fg given its early-post-i?-maximum 
spectrum is fairly normal, and so dominated by IME lines). 

Given two values of each of Pcentrai,WD (i-e., 2.9 x lO^gcm^^ and 10^°gcm~^), g (i.e., 
0.7 and 1), a (i.e., 1.2 and 1.7), and h (i.e., 1 and 0), there are 16 possibilities which we 
represent by the 8 SN 2003fg points on the Figures 1 and 2: N (^ = 1, a = 1.7), O (c/ = 1, 
a = 1.2), X (^ = 0.7, a = 1.7), and Y (^ = 0.7, a = 1.2). 

The input parameter values for all the SN 2003fg points and the lower bound masses 
we derive from them are given in Table 1. Table 1 also gives other output parameter values 
for the input parameter value sets: i.e., fco, /ime, /ipe, /ssni (which is the mass fraction of 
^^Ni), Ve, and E. The lower bound masses are actually obtained by solving SSC model using 
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MiPE-core, 'fiPE-core, Pcentrai,WD , and /i as 4 independent parameters and the procedure given 
in Appendix C. The last-hne parameter values with vipE-core = 8000 km s~^ give no solution 
for the mass since even a non-rotating WD for the given input parameters is too massive to 
yield this t'ipE-corc value. So we set fipE-corc = 7639.5 km s~^ which is the largest fipE-corc for 
which a solution exists given the other input parameter values. The output mass is just the 
mass of a non-rotating WD. The output masses of Table 1 can also be obtained, of course, 
by interpolation from Figures 1 and 2. 

Wc regard all the output masses in Table 1 as lower bounds on the actual mass since, as 
emphasized above, 8000 km s~^ is an upper bound on the IPE-core velocity and 7402.947 km s~ 
is cither an upper bound or at least not far removed from an upper bound. 

The lower bound masses in Table 1 are rated by the number of extreme parameter values 
(EPVs) used to calculate them: the lower the number, the more probable we think the lower 
bound mass. The most probable set of parameters is for point Y on Figure 1 (0 EPVs) and 
the least probable is for point N on Figure 2 (4 EPVs) for which in fact there is actually no 
M solution for viPE-core — 8000 km s"^ as noted above. The least probable set of parameters 
for a point on Figure 1 are for point N (3 EPVs). 

We see that the least probable lower bound masses for each figure (which are those 
associated with the N points) are approximately the Chandrasekhar mass. Excluding these 
lower bound masses gives a lowest lower bound mass of 1.62 M©. Excluding lower bound 
masses with two or more EPVs leaves 5 lower bound masses: the lowest of these is 2.16 Mq. 
The formally most probable lower bound mass (associated with point Y on Fig. 1) according 
to om rating system is the highest one, 2.84 Mq. This lower bound mass is close to the mass 
estimate of ~ SMq for SN 2003fg that we obtained from the simple-scaling argument (§ 2). 

From the above values we conclude that it is very probable that SN 2003fg is super- 
Chandrasekhar. The only lower bound masses that are approximately the Chandrasekhar 
mass are for extreme and improbable sets of the parameters. The more probable sets give 
lower bound masses that are super- Chandrasekhar and the 5 most probable sets give lower 
bound masses that are > 2 Mq. Given that YL's binding energy formula is only verified for 
M in the range ~ 1.4-2.05 Mq, a SN 2003fg mass > 2 Mq is as much as one can say based 
YL's binding energy formula for rotating super-Chandrasekhar-mass WDs (Appendix B) and 
the SSC model (§ 4). Given that the analysis of Langer et al. (2000) of accretion histories of 
WDs in binary systems limits rotating WD mass to < 2 Mq, our analysis leads to the final 
conclusion that the mass of SN 2003fg is probably ~ 2 Mq. 

The limitation on mass from Langer et al. (2000), has imphcations for the possible 
parameter values for g and a. In order to have the SN 2003fg mass ~ 2 Mq, the product of 
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ga would have to be more than the fiducial value of 0.7 x 1.2 = 0.84. Demanding M — 2 M©, 
solving for MipE-core for this M using the vipE-core = 8000 km and the Pcentrai.WD and h 
values of this section and Newton-Raphson method procedure of Appendix D, and using 
H2006's value for -£'boi/-E'56Ni(^boi) of 1.55 M© allows us to determine ga from 

ga — -. (24) 

-^IPE-core-E'SSNi (ibol) 

which we obtain using equation (6) of § 2. The ga values we obtain are 1.33 (pcentrai,WD = 

2.9 X lO^gCm-^, h = 1), 1.07 (Pcentral.WD = 2.9 X lO^gCm-^, h = 0), 1.21 (pcentral.WD = 

10^°gcm~^, h = 1), and 1.00 (pcentrai,WD = lO-'^'^gcm"^, h = 0). If we also demand that 
g = 0.7 (which as we argued in § 2 seems a good representative fiducial value), then the 
following a values are needed 1.91 (pccntraLWD = 2.9 x lO^gcm"'^, h = 1), 1.53 (pccntraLWD = 

2.9 X 109gCm-3, h = 0), 1.72 (pccntral,WD = lO^OgCm'S, h = 1), and 1.44 (pecntral,WD = 

10^°gcm~'^, h = 0). Since a is really not a well constrained parameter, these a values are 
not implausible. Only one is significantly outside of the semi-empirical a range of ~ 0.6-1.7 
we found in § 5. Recall that range is not definitive. Three of the a values are outside of the 
a range of 0.62-1.46 suggested by the SN la models of Hofiich & Khokhlov (1996). Recall 
that range was only suggestive. One must remember that SN 2003fg is a very peculiar and a 
unique SN la, and so its a value is not well constrained by the possible a values of existing 
SN la models and other known SNe la. 

One last fine point to discuss is related to the gamma-ray-escape effect discussed in 
§2. It follows from the discussion in § 2, that the gamma-ray-escape effect should give a 
tendency to lower a as /seNj increases if that increase brings ^^Ni matter closer to the surface 
as, in fact, it does in the SSC model. From Table 1, we see that the parameter sets with 
q; = 1.2 have larger /se^i values than the counterpart sets with a = 1.7 in all cases. Thus, 
coincidentally (since we are not invoking gamma-rays at all in the SSC model) a varies with 
/56Ni in a way that is qualitatively consistent with tendency of the gamma-ray-escape effect. 



7. CONCLUSIONS AND DISCUSSION 

Our analysis of the data for SN 2003fg supports the conclusion of II2006 that this SN la 
is supcr-Chandrasekhar (§ 6). We find it very improbable, although not absolutely ruled out, 
that SN 2003fg has only about a Chandrasekhar mass (the physical Chandrasekhar mass 
being 1.38 M© Appendix A) and probable that its mass is ~ 2 Mq as H2006 also concluded. 
Our conclusion relies on the adequacy of YL's binding energy formula for rotating super- 
Chandrasekhar-mass WDs (Appendix A), the analysis of Langer et al. (2000) that limits 
WD mass in possible accretion histories to < 2Mq, and our SSC model of SN la ejecta 



-28- 



(§4). It is likely that improvements to YL's binding energy formula will not change our 
SSC model results qualitatively, unless rotating WDs with masses significantly greater than 
the Chandrasekhar mass are ruled out by improvements in the input physics of YL's theory 
of differential rotation in mass- accreting WDs. 

An ultimate verification that SN 2003fg is super-Chandrasekhar would require that 

only realistic explosion models of super-Chandrasekhar-mass WDs allow a radiative transfer 
history (calculated with detailed radiative transfer) consistent with observations. Since we 
do not yet have completely adequate explosion models for any SN la, it will probably be 
some time before this ultimate verification is possible. 

What is the significance of SN 2003fg? Given that as of 2006 September 28 there are 
1347 SNe la or possible SNe la that have been observed according to the Central Bureau for 
Astronomical Telegrams (2006) (with, of course, a wide range of observational coverage from 
almost nothing to extensive) and SN 2003fg is the first SN la discovered for which there is 
significant evidence of a significant super-Chandrasekhar mass, we conclude that SN 2003fg- 
hke SNe la must be rather rare. Thus, for cosmological evolution and for cosmological 
distance determinations, SN 2003fg-hke SNe la are probably of small direct importance. 

But for the study of SNe la, SN 2003fg may be quite important. As YL point out, it 
is somewhat surprising that super-Chandrasekhar SNe la have not been much considered 
given that WDs rotate, and thus can be super-Chandrasekhar especially if their rotation 
is accelerated by accretion. Of course, the lack of observational need to consider super- 
Chandrasekhar SNe la has been a main factor in the relative neglect of such events. Also 
the relative homogeneity of SNe la pointed to a relative of homogeneity progenitors which 
pointed to maybe only near- Chandrasekhar mass SNe la. However, it could be (as YL and 
Tornambe & Piersanti (2005) have discussed at greater length) that there is a distribution of 
SN la progenitor masses from slightly below the Chandrasekhar mass to significantly above 
it and that the distribution of masses may explain part of the dispersion of SN la behavior. 
This hypothetical distribution is probably fast declining with mass coordinate in the vicinity 
of the probable mass of SN 2003fg given that SN 2003fg-hke SNe la seem so rare. In fact, 
it is possible that having a small super-Chandrasekhar mass is an important ingredient in 
normal SNe la explosions and radiative transfer histories. This is a speculation, but one that 
is hard to ignore in light of SN 2003fg. 
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0506028, and the Homer L. Dodge Department of Physics & Astronomy of the University 
of Oklahoma. We thank S.-C. Yoon for answers to questions about his work and Amedeo 
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A. 



THE CHANDRASEKHAR MASS 



The expression Chandrasekhar mass is used to mean the upper hmit on the mass of a 
non-rotating WD. Above this upper hmit, collapse occurs to a neutron star caused by the 
inability of the pressure of the matter's equation of state (EOS) to sustain the WD against its 
self-gravity. But what exactly the Chandrasekhar mass is depends on the degree of physical 
realism specified. 

Chandrasekhar in his book on stellar structure (Chandrasekhar 1957, p. 412ff) derives 
the upper limit based on Newtonian gravity and the zero-temperature, perfect Fermi electron 
gas in the extreme (special) relativistic limit where the EOS is 



where P is pressure, X is a constant depending on fundamental constants and the compo- 
sition of the WD, and p is (mass) density. It is the electron gas pressure that supports the 
WD against self-gravity. The WD with this EOS is an index-3 polytrope with a fixed density 
profile: polytropes being hydrostatic spheres with EOSs of the form 



where K is a. constant and n is the polytropic index (e.g., Chandrasekhar 1957, p. 84ff; 
Clayton 1983, p. 155fF). From Chandrasekhar's derivation (also given by, e.g., Shapiro 
& Teukolsky (1983, p. 61-64)) with the CODATA 2002 fundamental constant values (e.g., 
NIST 2006) and the modern solar mass value Mq = 1.9891 x lO^^g (e.g.. Cox 2000, p. 12), 
one obtains for the mass limit what can be called the formal Chandrasekhar mass. 



where the limitation to five significant digits for the coefficients is set by the number of 
significant digits given for the modern value of the gravitational constant {G — 6.6742(10) x 
10~^ in CGS units) and for the modern value of the solar mass. (We used double precision 
arithmetic to calculate the coefficients for equation (A3) and for all the calculations in this 
appendix.) The Chandrasekhar mass is independent of density. Thus, we can imagine scaling 
up the density throughout the WD by a constant factor and make the electron gas more 
and more relativistic throughout the WD, and thus make the extreme relativistic condition 
of derivation more and more exact. Thus, Chandrasekhar mass is asymptotically exact for 
zero-temperature, perfect Fermi electron gas as the density is scaled up toward infinity. The 
radius of the WD goes asymptotically to zero as Pcentrai,WD oo (e.g., Shapiro & Teukolsky 
1983, p. 63). Thus, if mass accretion increases a WD mass toward the Chandrasekhar mass. 



P = Kp'/' , 



(Al) 



(A2) 




(A3) 
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the WD in the original derivation picture is compacting toward a point mass. At some point 
in the accretion, one would expect the input physics to break down and some process would 
prevent compaction to a point mass. 

The lie quantity in equation (A3) is the mean molecular mass per electron and is given 



where the sum is over all elements i in the composition, Xj is the mass fraction of element i, 
Zi is the atomic number of clement i (we assume that the atoms are all completely ionized 
and the electrons all behave as free particles), and Ai is the atomic mass of element i in 
atomic mass units (AMUs) (e.g., Clayton 1983, p. 84). The effect of composition on the 
Chandrasekhar mass enters only through /Xg. For compositions not containing hydrogen or 
heavy elements Zi/Ai ~ 1/2, and so the fiducial value of //e is 2. For a CO WD with equal 
parts carbon and oxygen and standard atomic masses (which are some kind of terrestrial 
weighted averages of isotopic atomic masses) (e.g., NIST 2005), He — 2.00085 and Mch — 
1.4546 M0. 

In evaluating the formal Chandrasekhar mass, Chandrasekhar used older fundamental 
constant values and an older solar mass value, of course. He also used the proton mass where 
a modern person would probably use the AMU in his formulae, and this required that his 
atomic masses for evaluating /ig are in units of the proton mass rather than in units of the 
AMU (Chandrasekhar 1957, p. 415-416,432-433). For the first coefficient for equation (A3), 
Chandrasekhar obtained 5.75 (Chandrasekhar 1957, p. 423) from which the second coefficient 
value 1.44 follows. The coefficient values 5.75 and 1.44 appear in various places in the 
literature (e.g., Ostriker & Bodenheimer 1968). Using the fundamental constant values and 
solar mass value ostensibly used by Chandrasekhar (1957, p. 487) and Chandrasekhar 's 
definition of /le, we obtain for the coefficients 5.7513 and 1.4378 where we quote more digits 
than are physically significant to allow consistency checks: Chandrasekhar correctly rounded 
off to 3 digits given the precision of the values he used. If one does the calculations for the 
coefficients using the EOS parameter values calculated by Chandrasekhar (1957, p. 416) and 
his gravitational constant (G = 6.67 x 10^^ in CCS units: Chandrasekhar 1957, p. 487) rather 
than using the older fundamental constant values (except for the gravitational constant) 
directly (which implies using Chandrasekhar 's definition of He) and using the modern solar 
mass value, one obtains for the coefficients 5.7444 and 1.4361 where we quote more digits than 
are physically significant to allow consistency checks. If we used the modern fundamental 
constant values and solar mass value and Chandrasekhar 's definition of /le, the coefficients 
in equation (A3) would be 5.7395 and 1.4349 to the correct number of significant digits. We 
report the numbers in this paragraph to illustrate some of the sources of small inconsistencies 
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in formal Chandrasekhar mass values that may appear in the literature. 

The formal Chandrasekhar mass is an ideahzation in several respects. If one makes a 
correction for general relativity, then there is a collapse caused by gravity that occurs at a 
finite density which when considered apart from other effects is 2.646 x 10^°(//e/2)^ gcm~^ 
(e.g., Shapiro & Teukolsky 1983, p. 160) and at mass shghtly different from that of equa- 
tion (A3). Since the collapse occurs at a finite density, the electron gas is not in the extreme 
relativistic limit throughout the WD and gives rise to another correction for the collapse 
mass. The two effects reduce the mass of collapse from Mch,formai to 1.418 (2/yUe)^ Mq (e.g., 
Woosley & Weaver 1994, p. 84). There is also a Coulomb correction to account for the 
concentration of positive charge into nuclei rather spread uniformly in the electron gas (e.g., 
Shapiro & Teukolsky 1983, p. 29fF). Adding this effect gives a corrected Chandrasekhar mass 



(Shapiro & Teukolsky 1983, p. 32; Baron & Cooperstein 1990; Woosley & Weaver 1994, 
p. 84), where ctfs = 1/137.03599911(46) is the fine structure constant (e.g., NIST 2006) and 
where we have approximately generalized the expression from the references by replacing 
atomic number Z by mean atomic number (Z). For a CO WD with equal parts carbon and 
oxygen, //g = 2.00085 again, and {Z) — 7, one obtains Mch — 1.38 M©. 

Further corrections to the corrected Chandrasekhar mass are needed for finite temper- 
ature (Baron & Cooperstein 1990; Woosley & Weaver 1994, p. 85-86). These corrections 
increase the corrected Chandrasekhar mass, but for CO WDs they do so by only of order 



In fact, it is beheved that non-rotating CO WDs can never reach the exact corrected 
Chandrasekhar mass (whatever that exactly is) for a density-insensitive EOS. It is believed 
that WDs can only grow close to the corrected Chandrasekhar mass by accretion. During 
this process when central density Pcentrai.WD reaches of order lO^^gcm^^, electron capture on 
nuclei to make neutrons is likely to induce collapse to a neutron star by rapidly diminishing 
the electron gas pressure (Nomoto & Kondo 1991; YL). This collapse event depends on the 
central temperature being low enough, but that seems likely in likely accretion histories. If 
the central temperature is relatively high (of order 2.5 x 10^ K), then unstable carbon burning 
is hkely to lead to thermonuclear explosion when the central density is of order 3 x 10^ g cm"^ 
(e.g., Woosley & Weaver 1994, p. 108). The central density and temperature for collapse or 
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explosion actually do depend on the accretion history and are also are subject to revision 
with improved input physics. However, at present the non-rotating WD mass at the time 
of explosive carbon ignition is found to be about 1.38 (e.g., Woosley & Weaver 1994, 
p. 108; Nomoto et al. 1984): i.e., close to the zero-temperature corrected Chandrasekhar 
mass Mch. In the main text, we will just refer to 1.38 M© as the Chandrasekhar mass or 
physical Chandrasekhar mass for brevity. 

As a necessary simplification for their rotating WD modeling, YL adopted the zero- 
temperature, perfect Fermi electron gas EOS derived by Chandrasekhar in 1935 (Chan- 
drasekhar 1935). This EOS is sometimes called the Chandrasekhar EOS (e.g., Ostriker 
& Bodenheimer 1968). The Chandrasekhar EOS is more advanced than that used to derive 
the formal Chandrasekhar mass in that it allows for the full range of special relativistic effects 
from none to extreme. YL actually use the parameter values for the Chandrasekhar EOS 
that are given Chandrasekhar 's book (Chandrasekhar 1957, p. 416) and adopt /Xe = 2 for 
all calculations. By implication, YL should be using Chandrasekhar 's definition of fie- With 
this definition, fi^ should be 1.98640 for equal amounts of carbon and oxygen (calculated 
using standard atomic masses times the AMU divided by the proton mass), not 2. YL's 
value for the Chandrasekhar mass for their non-rotating model with infinite central density 
is 1.436 M0 which is just the Chandrasekhar mass value one would calculate using Chan- 
drasekhar's EOS parameter values, Chandrasekhar's value of the gravitational constant, the 
modern solar mass value, Chandrasekhar's definition of /le, and /le — 2. YL's Chandrasekhar 
mass value was, in fact, determined as just described (Yoon 2006). 

Because YL have omitted the general relativistic. Coulomb, and temperature correc- 
tions discussed above to the EOS, their masses will likely be in error by a few percent for 
the given input parameters. The Coulomb correction is probably the most important omis- 
sion. It tends to reduce mass for a given central density. Thus, YL's masses are probably 
overestimates for their specified central densities. Actually, YL's treatment of the EOS (us- 
ing older parameter values for the EOS and /Xg = 2 as described just above) may partially 
compensate for omitting the Coulomb correction. For example, the Coulomb correction re- 
duces the Chandrasekhar mass by about 2.5 % and and YL's treatment reduces the formal 
Chandrasekhar mass for equal amounts of carbon and oxygen from 1.4546 (see above) 
to 1.436 M0 or by about 1.3%. 
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B. THE YOON & LANGER BINDING ENERGY FORMULAE FOR 
ROTATING SUPER-CHANDRASEKHAR CO WHITE DWARFS AND 

SOME RELATED RESULTS 

The YL formula for the binding energy of rotating super-Chandrasekhar CO WDs is 

-E'bind(Pcentral,WD, M) — -Ebind,NR(Pcentral,WD) 

+-£'bind,coef(Pcentral,WD) [M — MNR(pcentral,WD)] ' (Bl) 

(YL's eq. (32)), where the energy units are foes (i.e., units of 10^-^ ergs), Pcentrai,WD is the 
WD central density in grams per centimeter cubed, M is the WD mass in solar mass units. 



^bind,NR(Pcentral,WD) = 0.1 X <{ -32.759747+6.7179802 X l0g(pcentral,WD) 

-0.28717609 X [log(pcentral,WD)]4 (B2) 



is the binding energy in foes for a non- rotating WD with central density Pcentrai,WD (YL's 
eq. (34)), 

^bind,coef(Pcentral,WD) = 0.1 X | -370.73052 + 132.97204 X log(pcentral,WD) 

-16.117031 X [l0g(peentral,WD)]' 

+0.66986678 X [log(pcentral,WD)]'| (B3) 

in foes (YL's eq. (33)), and 



-^NR(Pcentral,WD) = 1-436 X 



1 - exp<j -0.01316 X [l0g(pcentral,WD)]' '°' 
+0.2493 X log(pcentral ,WD 



(B4) 



is the mass in solar mass units of a non-rotating CO WD with central density Pcentrai,WD 
(YL's eq. (22)). The binding energy formula is a fit to the numerical results of YL and is 
verified for central densities in the range 10^-10^° gcm~^ and mass in the range ~ 1.16- 
2.1 Mq. Recall from § 2 that the CO WD central density range for explosion as SNe la 
is ~ 2.5 X 10*^-10^° gcm-^ (Woosley et al. 2006; Nomoto & Kondo 1991; YL). We have 
extrapolated the use of the binding energy formula up to 3Mq. 

The formula for MNR(pcentral,WD) is verified for Pcentral.WD > lO^gCm"^: for Pcentral,WD = 

2.5 X lO'^gcm-^, it gives 1.384 Mq and for Pcentrai.WD = 10^°gcm-^, 1.414 M©. Note that 
as Pcentrai.WD gocs to infinity, MNR(pcentrai,WD) gocs to 1.436 Mq. The value 1.436 Mq is the 



formal Chandrasekhar mass for YL's choice of equation of state and its parameters (see 
Appendix A). For the purposes of the investigations of this paper, the differences between 
1.436 Mq and the physical Chandrasekhar mass of 1.38 (see Appendix A) and between 
YL's equation of state and the exact equation of state are small. 

The binding energy formula is, in fact, almost linear in M for a fixed Pcentrai.WD (see 
YL's Fig. 9). Doing least-squares fits of the binding energy formula to a line for the mass 
range ~ 1 .4-3 Mq for the limiting central densities for SNe la give 

^bind(Pcentrai,WD = 2.5 X lO'' gcm-^ M) = (-1.133 ± 0.008) + (1.169 ± 0.003) X M (B5) 

and 

^bind(Pcentrai,WD = 10^° gcm-^ M) = (-1.923 ± 0.011) + (1.746 ± 0.005) X M , (B6) 

where again the energy is in foes and the mass in solar mass units. For exploring the behavior 
of the formula for kinetic energy of a SNe la (eq. (18) in § 4), it is, in fact, convenient to use 
a linear fit to the binding energy. With such a fit, the kinetic energy formula is 

E « ^ipeM[/ipe + Rh{l - /ipe)] - a - bM , (B7) 

where from the above we know that intercept a and slope b vary with allowed central densities 
for SNe la over -1.133 to -1.923 and 1.169 to 1.746, respectively. 

One use of the linear approximation to the binding energy formula is to verify that 
is roughly constant for fixed fractions of explosion-synthesized elements. We made 
use of this constancy in § 2 in our simple scaling argument for SN la mass. Substituting in 
EiPE = 1.61foe/M0 and R — 0.768 (see § 2) into equation (B7) and dividing that equation 
by M gives 

^ ^ 1.61 X [/ipe + 0.768 X h{l - /ipe)] - ^ + ]^ ■ (B8) 

To have any kind of explosion at all E/M > 0, and so the constant term of equation (B8) 
must be > — |amin|/M;h ~ —1.9/1.4 fa —1.4 (where amin = —1.923 is the minimum value of 

a for the Pccntrai.wD allowed range). If wc choose /ipe = 0.648, h = 0.798, and Pcentrai,WD = 
2.6 X lO^gcm^"^ (which are the values of model W7 (Nomoto ct al. 1984; Thielemann et al. 
1986)) as representative of SN la explosions, then a = —1.15 and b = 1.18 (from a least- 
squares fit to the binding energy formula over the mass range ~ 1.4— 3 iVfQ with Pccntrai,WD — 
2.6 X lO^gcm"^) and the first term of equation (B8) is 1.39. (Note /ipe = 0.648 and 
h — 0.798, imply fco = 0.0711 and /ime = 0.281 using equation (17) in § 4.) Substituting 
the calculated values for a, 6, and the first term of equation (B8) into equation (B8) gives 
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Thus, with the representative values, an explosion can be expected for any M since the 
constant term of equation (B9) is positive, but the value of E/M can only be very crudely 



and we can accept that it is so crudely. (Making use of equation (B9) and converting units, we 
find for M going from 1.38 M© to 3 Mq that ^J2EJM goes from 10260 km s"^ to 7750 km s'^) 
This constancy verifies one of our assumptions for the simple scaling argument provided that 
the SSC model is vahd up to ~ 3 Mq which requires that the usage of YL's binding energy 
formula can be extrapolated to ?>Mq. 

C. SOLVING THE SSC MODEL FOR MASS TREATED AS A 

DEPENDENT PARAMETER 

To calculate Figures 1 and 2, we used the SSC model mass M as an independent 
parameter. (See §§ 2 and 4 for a description of the parameters and the SSC model.) One 
can make M a dependent parameter and replace it among the independent parameters with 
■f^iPE-core- The solutiou for M can then be done via the Newton- Raphson method (e.g.. 
Press et al. 1992, p. 355ff) with IPE-core mass MipE-core, IPE-core velocity I'iPE-core (now 
independent), central density PccntraLWD, and h as the 4 independent parameter values. We 
assume the independent parameter values are given. The dependent parameters besides M 
are the e-folding velocity Vg, kinetic energy and the mass fractions /co, /ime, and /ipe- 

The Newton- Raphson method starts with an initial input mass Mi (for which we suggest 
an initial value below) and iterates to improve the input values toward the true value for 
M. The steps of the iteration are as follows. An input Mj for the ith iteration along with 
the given AfipE-corcj Pccntrai,WD and h values are used to evaluate the following version of the 
kinetic energy formula (see eq. (18) in § 4): 



where we have used the fact that MipE-corc = ^/ipe (sec cq. (4) in § 2) to eliminate /ipe 
from equation (18). Using the calculated ith iteration kinetic energy E'j, we obtain the ith 
iteration e-folding velocity 



where we have used equation (12) (§ 4), but with assumption that E^ is in foes and Mj is in 
solar mass units, and so need conversion factor B — lO^-*" ergsfoe^-^/MQ to give v^^i in CCS 
units. The ith iteration IPE-core reduced velocity xiPE-core,* = '^iPE-core/'i^e.i (see eq. (14) in 
§ 4) (assuming the given t'lPE-core value is in CCS units). One then uses equation (15) (§ 4), 




E — £'lPE[-MiPE-core + Rh{M — MiPE-core)] " -Sbind (Pcentral,WD , M) , 



(CI) 




(C2) 
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to find the ith iteration /iPE,i and then the ith iteration output mass Mout.i of the iteration 
follows from 

MipE -core 



out,j 



IPE,i 



(C3) 



where wc have made use of equation (4) in (§ 2) For the Newton-Raphson method iteration, 
the quantity whose zero we want to find as a function of M is M^ut — M. The ith iteration 
Newton-Raphson correction AMj to approximately zero M^ut — M is given by 



Mout - M 



dMout/dM - 1 



(C4) 



Mi 



It is straightforward to show that the ith iteration value of the derivative dMo^^/dM is given 

by 



dM, 



out 



dM 



Mi 



MipE-core df_dx_dVe dw_ \ 

dx dVf, dw dM I 



P 

MiPE- core 



Mi 



IPE,i 



X 



IPE-core,i ' 
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3^IPE-core,i) 




'd{E/M)' 










dM 




Mi 



, (C5) 



where w = {B/6){E/M), 



d{E/M) EiPERh - dEund/dM 



dM 



M 



E 



and 



dE^ 



bind 



dM 



1.03 X £'bind,coef(Pcentral,WD) [M — MNR(pcentral,WD)] 



0.03 



(C6) 



(C7) 



The last equation was obtained by differentiation from equation (Bl). 



The Newton-Raphson method in this case converges quite quickly in practice. The 
guarantee of convergence is that Afout ~ M is strictly decreasing with respect to M with 
no stationary points for the allowed Pccntrai.wo) range for SN la explosion of ~ 2.5 x 10^- 
lO^^gcm^^ (Woosley et al. 2006; Nomoto & Kondo 1991; YL). This means that the Newton- 
Raphson correction AMj always has the correct sign and if the corrections are prevented 
from causing any overshooting of the diminishing allowed range for the true mass M, then 
convergence must follow. All that is required to show the strictly decreasing nature of 
Mout — M is to show that E/M is monotonically decreasing with respect to M with a 
stationary point only at M = oo. This is clear from equation (C5) since all the factors on 
the right-hand side of the second expression are clearly always positive or zero except for 
d{E/M)/dM. (Note that Mout - M has no stationary points even though d{E/M)/dM and 
dMout/dM do because the — M term derivative is —1.) 
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Because of the unusual exponent in the binding energy formula (see eq. (Bl) in Ap- 
pendix B), it is awkward and not enlightening to work with the exact formula for kinetic 
energy equation (CI). Instead we will show that E/Al is monotonically decreasing with 
respect to M using the kinetic energy formula with the linear binding energy formula ap- 
proximation: sec equation (B7) in Appendix B. First, we rewrite equation (B7) in terms of 
the proper variables for the Newton-Raphson method: 

E ^ £;iPE[MiPE-core + i?M^-MPE-core)] 

= £;iPEMipE-core(l - - O + (£;iPEi?/i - &)M , (C8) 

where we note that the intercept term £'ipEMipE-core(l — a > always since 1 — Rh > 

always and a < for all linear fits to the binding energy formula for its allowed Pcentrai.wo) 
range (see Appendix B). The sign of the slope Ei^^Rh — b depends on the relative sizes of h 
and b. Using the model W7 (Nomoto et al. 1984; Thielemann et al. 1986) values h = 0.798 
and b = 1.18, we find that Ei^^Rh — b — —0.19: thus, what we can take as the representative 
slope of equation (C8) is negative. Using equation (C8), we obtain 

I ^ E,,M.^^Ul -Rh)-a ^ ^^^^^^ _ ^ ^^g^ 
and then differentiating E/M we get 



d{E/M) 
dM 



■glPEMiPE-core(l " Rh) - O 
M2 



(CIO) 



Since £'iPE-/^iPE-core(l — Rh) — a > always, it is clear within our approximation for E that 
E/M is monotonically decreasing with the only stationary point at M = oo. Now wc have 
shown this behavior only for the approximate equation (C9). But since this equation is, in 
fact, a very good approximation to equation (CI) divided by M and the behavior is robust for 
all allowed values of a (sec Appendix B), we conclude that £'/M is to at least good accuracy 
monotonically decreasing with M in general with a stationary point only aX M = oo. It then 
follows from equation (C5) that Mout — M is to at least good accuracy strictly decreasing 
with M. This is what we needed to show to prove that the Newton-Raphson method would 
always converge for the allowed Pcentrai.wo) range. 

We can deduce some approximate limiting behaviors of the Mout function by making 
use of equation (C9). If EipeRH — b < (which we take to be the representative case) and 
M is increased, then at some point (which is at finite M if Eip^Rh — 6 < and infinite M 
otherwise) E/M — > 0, — > 0, x — > oo, / — > 1, and Mout MipE-core- On the other hand, 
if Ejp^Rh — b > and M — > oo, then E/M, Ve, x, /, and Mout approach asymptotically to 
a finite positive value, a finite positive value, a finite positive value, a finite positive value 
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less than 1, and a finite positive value greater than MiPE-core, respectively. If M — > 0, then 
E jM ^ oo, Ve — > oo, a; ^ 0, / — > 0, and Mo^t oo. Actually, we cannot let M — > in the 
Mout function since YL's binding energy formula has no solution for M < Mnr. 

What is the appropriate initial value Mi for the Newton-Raphson method? Well the 
lower bound on the true value for M is clearly 



since a physical M cannot be less than the input core mass MipE-corc nor than the mass 
of a non-rotating WD for the input density MNR(pcentrai,WD)- In practice, we have found 
-Slower bound to work wcU as an initial value Mi. 

If we use Miower bound aS an input mass M and find that Mout (Mower bound) < Mower bound, 

then there is no solution since Newton-Raphson correction from equation (C4) will be neg- 
ative. Recall that Newton-Raphson correction always has the correct sign for our case, and 
so a negative correction to the lower bound on M implies there is no solution. One way 
of looking at this case is to say there is not enough nuclear burning energy released with 
any allowed value of M to match the input t'lPE-core- On Figures 1 and 2, the no-solution 

case with Miower bound = MNR(pcentral,WD) and Mout (Miower bound) < MNR(pcentral,WD) WOuld 

correspond to having the point (MiPE-corc, "^^iPE-core) in the region to the left of the Mnr mass 
curves: if was to the right, there would be a solution. This no-solution case is actually the 
case for the point N on Figure 2 for /i = (which means /ime = recall) (see § 6). 

If Miower bound givcs E < 0, then there is also no solution mass M. Since the intercept 
term of equation (C8) is always positive, then < for any input M implies the slope is 
negative. This means that the approximate kinetic energy function (and the exact kinetic 
energy function to high accuracy) can only grow more negative for M increasing beyond 
Miower bound- Thus for no allowcd M is the kinetic energy positive which is required for an 
explosion. There is not enough nuclear burning energy for any allowed M to explode the 
WD with the input parameter values in this case. 

The true values for /ire, v^, kinetic energy E, and mass M are obtained from the 
Newton-Raphson iteration. The true values of fco and /ime are obtained using the h and 
/iPE values in equation (17) (§ 4). 

Note that if we had used either of fco or /ime as an independent parameter instead of 
h, we would have to deal with an extra check on whether or not a given converged solution 
M was physically allowed. If we used fco (/ime) as an independent parameter and solved for 
M and then found that /jme < {fco < 0), the found solution M would be ruled out: i.e., 
there would be no physically allowed solution for the given input parameter values. Using 
h as an independent parameter eliminates this complication for valid h (i.e., h in the range 
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0-1) and a solution mass M > Miower bound, since one will only obtain valid fco and /ime 
values from equation (17) (§ 4): i.e., values in the range 0-1. Thus, the solution mass M 
can be known to be valid without having to check on the fco or /ime values. The above 
discussion shows why we prefer using h as an independent parameter to either of fco or 
/ime- 



D. SOLVING THE SSC MODEL FOR CORE MASS TREATED AS A 

DEPENDENT PARAMETER 

In this appendix, we show how to solve the SSC model for core mass MjpE-core treated 
as a dependent parameter. The independent parameters arc now mass M. IPE-core velocity 
■?7iPE-corc, WD central density Pccntrai,WD, and h. Wc assume the independent parameter 
values arc given. The dependent parameter besides MiPE-core are the e-folding velocity Ve, 
kinetic energy E, and mass fractions fco, /ime, and /ipe- The solution is by the Newton- 
Raphson method (e.g.. Press et al. 1992, p. 355ff) and is similar to that of Appendix C, but 
is somewhat simpler since the binding energy formula (eq. (Bl) in Appendix B) does not 
depend on MiPE-core- 

The Newton-Raphson method iteration starts with an initial input core mass MiPE-core,i 
(for which we suggest an initial value below) and iterates to improve the input values toward 
the true value for MipE-core- The steps of the iteration are as follows. An input MipE-core,? 
for the ith iteration is used to evaluate the ith iteration kinetic energy and e-folding velocity 
Ve,i using, respectively, equations (CI) and (C2). The ith iteration IPE-core reduced velocity 
a;iPE-core,i = WiPE-corc/^^e,j (scc cq. (14) in § 4) (assumiug the given ViPE-corc value is in CCS 
units). One then uses equation (15) (§ 4), to find the ith iteration /iPE,i and then the ith 
iteration output core mass MiPE-core,out,i follows from 

-^IPE-core,out,i = M fiPE,i , (Dl) 

where wc have used equation (4) (§ 2). For the Newton-Raphson method iteration, the 
quantity whose zero we want to find as a function of MiPE-corc is MiPE-core,out — AfiPE-corc- The 
ith iteration Newton-Raphson correction AMipE-core,i to approximately zero MipE-core,out ~ 
MiPE-core IS giveu by 



AMi 



IPE-core,i — 



MiPE -core, out MiPE -core 



(D2) 



dMipE-core,out/ dMip^- core ^ . 

It is straightforward to show that the ith iteration derivative dMiPE-core,out/c^-^iPE-core is 
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given by 



IPE-corc,out 



IPE-core 



M 



df dx dve dw 
dx dve dw dMiPE- 



^IPE-core,i ^^Pl 2;iPE-core,i j 
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X 



B 



EipE(l-m) , 



(D3) 



where w = {B/6){E/M). The derivative (iMiPE-core,out /'i^iPE-corc < always with the equal- 



ity holding only for xiPE-core,i = and xiPE-core.i = oo. (Note x[ 



3 

IPE-core 



,il e,i 



X 



IPE-core,i/ "^IPE-core 



Thus, the Newton-Raphson method is guaranteed to converge to the true MipE-core value 
since MiPE-core,out — -/^IPE-core IS strictly decreasing with MiPE-core provided the Newton- 
Raphson corrections are prevented from causing any overshooting of the diminishing allowed 
range for the converged solution and provided the independent parameter values allow a 
valid MipE-corc solution (i.e., one in the range 0-M). (Note that MipE-corc,out — ^iPE-corc has 
no stationary points even though MipE-core,out does because the — MipE-core term derivative 
is -1.) 

The true /ipe, e-folding velocity fe, kinetic energy and mass MiPE-corc values are 
obtained from the iteration. The true values of /co and /ime are obtained using the h and 
/ipe values in equation (17) (§ 4). 

A good initial input MiPE-corc,i for Newton Raphson iteration is M/2 which is just the 
midpoint of the allowed range 0-M. 
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Fig. 1. — Curves of constant SN la mass from the SSC model for super-Chandrasekhar-mass 
SNe la on a plot of IPE-core velocity versus IPE-core mass. There are two families of curves: 
the sohd curves are for only IMEs above the IPE core {h — 1) and the dotted curves are 
for only CO above the IPE core {h — 0). This plot is calculated for WD progenitors with 
central density 2.9 x lO^gcm"^. The plot and the data points are explained at length in the 
text (§§ 4, 5, and 6). 
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Fig. 2. — The same as Figure 1, except that the plot is calculated for WD progenitors with 
central density 10^°gcm~^. 



Table 1. LOWER BOUND MASSES AND OTHER OUTPUT PARAMETERS FOR SN 2003fg FOR VARIOUS 

CHOICES OF INPUT PARAMETERS FOR THE SSC MODEL 



Symbol'* 


MipB-core'' 


Pcentral.WD*^ 




a" 




No. ofs 




/go' 


/iME-" 


/ipb'' 


/56Ni' 




kinetic energy i?" 




{Mq) 


(108gcm-3) 








EPVs 












(kms-l) 


(foes) 


Y 


1.84 


2.9 


0.7 


1.2 


1.0 





2.84 


0.000 


0.350 


0.650 


0.455 


2390 


1.933 


Y 


1.84 


2.9 


0.7 


1.2 


0.0 


1 


2.44 


0.245 


0.000 


0.755 


0.528 


2024 


1.194 


X 


1.30 


2.9 


0.7 


1.7 


1.0 


1 


2.17 


0.000 


0.402 


0.598 


0.419 


2583 


1.731 


A. 


l-OU 


9 Q 


U. 1 


1 7 
1. ( 


U-U 


9 


1 

1 .o^ 




u.uuu 


U. / UD 




Ziyi 


1 .uoo 


O 


1.29 


2.9 


1.0 


1.2 


1.0 


1 


2.16 


0.000 


0.403 


0.597 


0.597 


2587 


1.726 


O 


1.29 


2.9 


1.0 


1.2 


0.0 


2 


1.83 


0.295 


0.000 


0.705 


0.705 


2195 


1.052 


N 


0.91 


2.9 


1.0 


1.7 


1.0 


2 


1.68 


0.000 


0.457 


0.543 


0.543 


2803 


1.571 


N 


0.91 


2.9 


1.0 


1.7 


0.0 


3 


1.40 


0.348 


0.000 


0.652 


0.652 


2383 


0.947 


Y 


1.84 


10.0 


0.7 


1.2 


1.0 


1 


2.52 


0.000 


0.269 


0.731 


0.512 


2104 


1.332 


Y 


1.84 


10.0 


0.7 


1.2 


0.0 


2 


2.28 


0.191 


0.000 


0.809 


0.566 


1840 


0.920 


X 


1.30 


10.0 


0.7 


1.7 


1.0 


2 


2.02 


0.000 


0.354 


0.646 


0.452 


2406 


1.392 


X 


1.30 


10.0 


0.7 


1.7 


0.0 


3 


1.77 


0.266 


0.000 


0.734 


0.514 


2095 


0.929 


O 


1.29 


10.0 


1.0 


1.2 


1.0 


2 


2.00 


0.000 


0.356 


0.644 


0.644 


2414 


1.394 


O 


1.29 


10.0 


1.0 


1.2 


0.0 


3 


1.76 


0.268 


0.000 


0.732 


0.732 


2102 


0.929 


N 


0.91 


10.0 


1.0 


1.7 


1.0 


3 


1.62 


0.000 


0.438 


0.562 


0.562 


2726 


1.437 


N 


0.91 


10.0 


1.0 


1.7 


0.0 


4 


1.414° 


0.356 


0.000 


0.644 


0.644 


2304 


0.896 



"Symbol is the SA'iiibol used for the data point on Figures 1 and 2. 

''This is the IPE-core mass. It depends inversely on the g and a parameters. In fact Mip^^^^^ga = 1.55 M©. We make the approximation that 
the IPE core contains all explosion-synthesized IPEs and that no explosion-synthesized IPEs are outside the IPE core. 

°The Pcentral.WD quantity is the central density of the progenitor CO WD. 

''The g fa<;tor is the ratio of ^®Ni mass to IPE mass which we also take to the IPE-core mass MiPB-core- 

"The a quantity is the bolometric-luminosity-maximum-light ratio of luminosity to the total instantaneous radioactive decay energy release rate 
per unit ^''Ni mass. 

^The h = /ime/(/co + /ime) quantity is the fraction of non-IPE element matter in the ejecta that is explosively-synthesized IME matter. 
^This is the number of extreme parameter values (EPVs) used in determining the lower bound mass. 
''The Mlb quantities are the lower bound masses for the given input parameters to the SSC model. 
'The fco quantity is the mass fraction of CO in the ejecta. 



j The /iMB quantity is the mass fraction of IMEs in the ejecta. 

'^The /ipE quantity is the mass fraction of IPEs in the ejecta and also interior mass fraction of the IPE core. 

'The /sej^ji quantity is the mass fraction of ''^Ni in the ejecta. It equals /ipe when g = 1. 

™Thc Ve quantity is the e-folding velocity of the ejecta. 

"This is the kinetic energy of the ejecta in foes (1 foe = 10^^ ergs). 

°The last-line input parameters give no solution for the mass for ■yiPE-core = 8000 km s~^ since even a non-rotating WD for the given input 

parameters is too massive to match this I'ipE-corc value. So we set iiipE-corc = 7639.5 km which is the largest ^'IPE-core for which a solution exists. 
The output mass is just the mass of a non-rotating WD. It is an upper bound mass or close to an upper bound mass. 

Note. — Sec the text § 2 for a full description the parameters (except h and /sej^;), § 4 for a description of the SSC model and the h parameter, 
and § 6 for a discussion of the parameter values and the table. The SSC model input parameters for the calculations of the lower bound masses are 
IPE-corc mass MipE-corei IPE-core velocity (8000 kms^^ for all cases, except for the last line where «ipE-core = 7639. Skms^-' was used), Pcentral WDi 
and h. The IPE-core mass was itself calculated from equation (6) (§ 2) using the g and h values and H2006's value of 1.55 Mq for i/bol/^56j,ji(tboi)- 
The quantities Mlbi /cOi /imEi /ipEi /seN;, Ve, and E axe output parameters of the SSC model. The output parameters were determined from the 
SSC model by the procedure given in Appendix C. The IPE-core velocity 8000 km s~^ is an upper bound on the actual IPE-core velocity, and so the 
output masses for this IPE-core velocity are lower bounds. 



